General Disclaimer 


One or more of the Following Statements may affect this Document 


• This document has been reproduced from the best copy furnished by the 
organizational source. It is being released in the interest of making available as 
much information as possible. 


• This document may contain data, which exceeds the sheet parameters. It was 
furnished in this condition by the organizational source and is the best copy 
available. 


• This document may contain tone-on-tone or color graphs, charts and/or pictures, 
which have been reproduced in black and white. 


• This document is paginated as submitted by the original source. 


• Portions of this document are not fully legible due to the historical nature of some 
of the material. However, it is the best reproduction available from the original 
submission. 


Produced by the NASA Center for Aerospace Information (CASI) 




^ s tems, ^ 


f NASA-CB- 135414-App) HELICCENTBIC 
INTEBPLANET ABY LOW THBUST TBAJECTOB Y 
OPTIMIZATION PECGBAM, SUPPLEMENT 1, FAST 2 
Final Bepott (Easiness and Technological 
Systems, Inc.) 103 p HC A06/MF A0 1 CSCI 22A G3/12 


N78-25 106 


Onclas 

2ieoi 


BTS-TR- 7 8- 5 4 


NASA CR-135414 


HILTOP SUPPLEMENT 

HELIOCENTRIC INTERPLANETARY LOW THRUST 
TRAJECTORY OPTIMIZATION PROGRAM 

SUPPLEMENT #1 


F.I. Mann 
J.L. Horsewood 


FINAL REPORT 

Part 2 of 2 
Contract NAS 3-20950 


February 1978 


prepared for 


NATIONAL AERONAUTICS AND SPACE ADMINISTRATION 
Lewis Research Center 
Cleveland, Ohio 44135 

by 


BUSINESS AND TECHNOLOGICAL SYSTEMS, INC. 
10210 Greenbelt Road, Suite 440 
Seabrook, Maryland 20S01 


SUMMARY 


This report is the first supplement to the currently existing 
primary HILTOP program document [11 and describes the improvements made 
to the HILTOP electric propulsion trajectory optimization computer 
program since the publication of the primary document. 

A new, more realistic propulsion system model has been implemented 
in the program, in which various thrust subsystem efficiencies and 
specific impulse are modeled as variable functions of power available to 
the propulsion system. The number of operating thrusters are staged, 
and the beam voltage is selected from a set of five (or less) constant 
voltages, based upon the application of variational calculus. The 
constant beam voltages may be -optimized individually or collectively. 

The new propulsion system logic is activated by a single program 
input key in such a manner as to preserve the old HILTOP logic. 

The report contains the new analysis describing these features, 
a complete description of program input quantities, and sample cases 
of computer output illustrating the new program capabilities. 
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NOMENCLATURE 


Generally:, upper-case symbols denote vectors and lower-case symbols 
denote scalars. Lower-case symbols with bars denote unit vectors. The 
abbreviations EPS for electric propulsion system and BVP for boundary value 
problem are used. 


a EPS instantaneous thrust acceleration, expressions (22) and (35); 

semi -major axis 




Semi -major axis of primary- target capture orbit 
Solar power law coefficients 

Arbitrary unit vectors used in (132) and (139) of [1] 

A coefficient in the efficiency law of the old spacecraft model; 
auxiliary quantity defined by (29) 

Launch vehicle coefficients 

Vector constant of optimal rocket problem, expression (63) of [1] 


C° Pvadians-to-degrees conversion factor 


Zeroth-order term in expansion for thrust reduction factor, 
expression (12) 

Cg Coefficient in first-order term in expansion for thrust reduction 

factor, expression (12) 


vii 


me blank 


Constant in expression (13) or (14) for n u 

Constant in expression (13) c-r (14) for n u 

Ratio of accelerator current to beam current, as used in 
(15) and (16) 

Zeroth-order term in expansion foi screen supply efficiency, 
expression (17) 

Coefficient in first-order term in expansion for screen 
supply efficiency, expression (17) 

EPS jet exhaust speed, constant in the "old” spacecraft model; 
defined by (20) or (36) or (56) in the improved model; 
abbreviation for cosine function 

Acceleration conversion factor, expression (35) 

Retro stage jet exhaust speed 

Speed conversion factor, expression (36) 

Auxiliary quantity given by expression (74) of [1] 

Coefficients in quadratic expression for Av.., expression (78) of [1] 
Digital Interface Unit 

A coefficient in the efficiency law; an auxiliary quantity in 
the coast-phase solution; solar flux density 

Eccentric arfomcily (a scalar) 


e A coefficient in the efficiency law; the base of the natural 

logarithms; eccentricity; subscript denoting Earth 

Spacecraft unit angular momentum vector 

"e Spacecraft unit radius vector 

EPS unit thrust vector 

e" Spacecraft unit velocity vector 

e x Retro stage characteristic speed exponential factor given by 

expression (76) of [11 

e^ Unit primer vector 

F Auxiliary scalar function defined by (215) of [1] 

F t See «F t 

F y Auxiliary quantity defined by (30) 

f EPS instantaneous thrust magnitude; f-f unction of the f and g 

series; subscript denoting a desired value; true anomaly; 
auxiliary variable defined by equation (147) of [1] 

f Retro stage thrust magnitude 

fj. Total thrust, expression (18) 

f Auxiliary quantity given by expression (77) of [1] 

A 

G Auxiliary quantity defined by (28) 

Auxiliary scalar functions in the coast-phase solution, equation 

(45) of [11 

ix 


g EPS reference thrust acceleration (old spacecraft model); 

g-function of the f and g series; BVP point-constraint 
geometric mean of the weighting factors; parameter local to 
expression (47) 

g 1 Parameter local to expression (47) 

g r Reference gravity acceleration constant, used in (19) 

g x Auxiliary quantity given by expression (97) of [1] 

H Spacecraft angular momentum vector 

h Magnitude of spacecraft angular momentum vector 

h Spacecraft unit angular momentum vector 

hj Component of h y containing thrust acceleration, expression (41) 

h y Variational hamiltonian, given by (54) 

h x | 

f Cartesian components of spacecraft angular momentum vector 
h z ) 

h^ Thrust-switching step-function (old spacecraft model only) 

Ig Beam current per thruster 

Igp Phantom beam current (per thruster), equal to the beam current 

a thruster would have if there were no limitations on Ig , 
given by (26) 

Ig^ Restricted phantom beam current, expression (35) 

I max Maximum beam current (per thruster) 


x 


Minimum beam current (per thruster) 


EPS specific impulse, expression (19) 

Reference value of I^p , defined in discussion following 
expression (23) 

Subscript pertaining to an intermediate target; inclination to 
ecliptic; general subscript or running index; inclination of 
parking orbit about Earth 

Unit vector along x-axis 

Parking orbit inclination associated with range safety limit 
Index-set of the BVP dependent variables 
Unit vector along y-axis 

Jettison indicator for solar arrays (or other power source) 
prior to retro maneuver, used in (5) 

Unspecified-reference-power indicator (old spacecraft model only) 

EPS propulsion system jettison indicator (retro maneuver) in old 
spacecraft model 

Retro stage existence indicator 

EPS tankage jettison indicator (retro maneuver) 

EPS thrust subsystem jettison indicator (retro maneuver), used in (5) 

Fundamental constant associated with Mercury propellent, used in 
(18) and (19); in [1], arbitrary positive constant associated with 
performance index (replaced by in this document); temporary 
variable ultimately equated to inverse of the characteristic 
degradation time 


k 



m 

in 


Unit vector along z-axis 

Auxiliary quantity given by expression (75) of [11 

Intermediate-target drop-mass factor defined by expression (6) 
of [1] 

Retro stage tankage mass factor defined by expression (11) of [1] 

EPS structure mass factor defined by expression (3) of [1] 

Intermediate-target sample-mass factor defined by expression 
(6) of [1] 

EPS tankage mass factor defined by expression (7) of [1] 

Launch site latitude (scalar) 

Mean anomaly (scalar) 


Coefficients used in computing nuclear and total magnitudes of a 
celestial body (scalars) 

Nuclear magnitude (scalar) 

Total magnitude (scalar) 

Spacecraft total mass variable 
Auxiliary unit vector given by expression (53) 
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m Mass flow rate, expression (21) 

m a Solar array (or other power source) mass, expression (2) 

Constant mass component of m g , expression (2) 
m drop Intermediate-target drop-mass given by expression (6) of [1] 
m £ Thruster neutral propellant loss, used in (13) and (14) 

m n Neutralizer propellant flow rate, used in (13) and (14) 

m net Net s P acecra f' t mass 

m Q Initial spacecraft mass (payload of launch vehicle) given by 

expression (2) of [1] 

nip EPS propellant mass 

m Electric propulsion system mass given by expression (4) of [1] 

m p Retro stage mass, expression (4) 

m rp Retro stage propellant mass given by expression (5), or (9) of [1] 

m rs Retro stage structure mass defined by expression (11) of [1] 

m rst Retro stage structure and tankage mass given by expression (11) of 

m 

m g EPS structure mass 

m samp Intermediate-target sample-mass given by expression (6) of [1] 
m.j. EPS tankage mass 


xiii 


m tti EPS thruster subsystem mass, expression (3) 

Anip Propellant mass increment due to primry- target spiral maneuver 

n Exponent in step-size law, expression (39) of [13; subscript denoting 

time at the primary target; number of BVP dependent variables 

n Unit vector normal to the solar arrays 

n max Maximum number of operating thrusters 

n min Minimum number of operating thrusters (>0) 

rip Unit vector directed along a planet's north pole 

n t Number of operating thrusters (variable) 

An^ Increment in n t , when staging thrusters 

o Subscript denoting launch time; subscript denoting the beginning 

of a computation step 

P A celestial body's position vector; BVP partial derivative matrix 

PDS Power Distribution System 

PPU Power Processing Unit 

p EPS instantaneous power in old spacecraft model; subscript 

denoting a perturbed, or neighboring, parameter; auxiliary 
variable in equations (79) of [1] 

Ap Ratio of housekeeping to reference power, Pf/P re f > old model only 

p Q Total instantaneous power developed by arrays (or other power source), 

expression (6) 


xi v 



p 

xonv 

p d 

p diu 

p h 

P^v 

p max 

p mni 

p ref 

p t 

Pts 

p to 

Pi 

p 2 


Maximum power output of power source that can be utilized by 
the thruster subsystem and other spacecraft modules, expression 
(25) 

Reference power of solar arrays or other power source, used in 
(2) and (6) 

Beam power, expressions (10) and (23) 

Housekeeping power in improved spacecraft model, expression (9) 
Power output of power distribution system, expression (7) 

Digital interface unit power requirement, used in (9) 

Housekeeping power in old spacecraft model 

Power processor low voltage power requirement, used in (9) 

Maximum power input to an individual thruster, expression (25a) 
Mission module power requirement, used in (9) 

EPS reference power (old model only) 

Power input each PPU, expressions (8) and (27) 

Thrust subsystem power requirement, used in (9) 

Reference power of each thruster, expressions (3) and (24) 

Auxiliary quantities in coast-phase solution, expressions (54) and 
(55) of [1] 
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s 


t 

*b 

Vatlo 

At 


U 


Auxiliary variable in equations (79) of [1]; solar array 
radiation damage factor 

Spacecraft position vector 

Magnitude of R 

Primary- target capture-orbit apocenter distance 

Earth-to-spacecraft communication distance 

Unit vector along line of ascending node 

Primary-target capture-orbit peri center distance; primary- 
target swingby passage-distance 

Swingby passage-distance unit vector 

Value of r for which y-curve is at a maximum 

Abbreviation for sine function; auxiliary variable used in 
equations (79) of [1] ; degradation time 

Unit vector directed toward Canopus 

Time 

Retro maneuver burn time given by expression (12) of [1] 
Minimum throttling ratio, expression (25b) 

Time-increment due to primary-target spiral maneuver 
Generalized universal anomaly during thrust phases 


Generalized universal anomaly increment, equivalent to the 
computation step-size during numerical integration 

Neutralizer to beam coupling potential, used in (15) and (16) 

i, L 

Beam voltage (I tn value selected from set of up to five constant 
values) 

Largest beam voltage 
Smallest beam voltage 

Discharge voltage, expressions (15) and (16) 

Hyperbolic excess velocity (or encounter velocity) 

Swingby planet arrival hyperbolic excess velocity 
Swingby planet departure hyperbolic excess velocity 
Magnitude of spacecraft velocity 
Characteristic speed of a rocket maneuver 
Escape speed from launch parking orbit 

Minimum velocity impulse required for non-coplanar injection from 
a circular orbit to a given excess velocity 

Speed of a spacecraft in a circular orbit 

Planetocentric speed at primary-target swingby closest-approach 
point; auxiliary speed given by equation (72) of El 3 
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Hyperbolic excess speed (or encounter speed) 


Retro stage imoulsive velocity increment magnitude; characteristic 
velocity associated with primary-target spiral maneuver; incremen- 
tal speed required at powered swingby 

Retro stage total velocity increment magnitude 

Minimum incremental velocity (magnitude) for coplanar boost out 
of circular orbit 

Velocity penalty due to noncoplanar boost out of circular orbit 

Velocity penalty due to launch azimuth 

Auxiliary variable in equations (79) 0 f [1] 

First Cartesian component of position; a general variable; a 
general state variable; auxiliary variable in equations (79) 0 f p] 

Second Cartesian component of position; auxiliary variable in 
equations (79) of [1] 

Third Cartesian component of position 

EPS specific mass; geocentric right ascension of launch excess 
velocity 

Thrust reduction factor due to double ions and beam divergence, 
given by (12) 


Auxiliary parameters defined by equations (211) and (212) of [1] 
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Specific mass of the solar arrays, used in (2) 


Communication angle (Sun-Earth-spacecraft) 

Specific mass of the power conditioning and thruster subsystem, 
used in (3) 


Arbitrary, independent angles defining orientation of excess 
velocity in (132) and (139) of [1] 

Independent variable of coast-phase solution, also generalized 
to be the independent variable on the entire trajectory 

Value of B at the beginning of a computation step 

Computation step size (increment of trajectory independent 
variable) 

Normalized power function 
ay/sr 

ay/9d, where d is the solar flux density 

Launch hyperbolic-excess-velocity asymptote declination; BVP 
dependent-variable tolerance 


Bend angles of hyperbolic arrival and departure trajectories, 
expression (213) of [1] 

Total bend angle given by expression (214) of [1] 

Kronecker delta function ORIGINAL PAGJ3 
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xix 


e Auxiliary quantity in the coast-phase solution; obliquity of 

the Earth's equator to the ecliptic 

Discharge losses, used in (15) and (16) 

n EPS efficiency (old model) 

n 1 dn/dc 

n conv DC- DC converter efficiency, used in (9) 

Hq Discharge supply efficiency, used in (16) 

Efficiency of the power distribution system, used in (7) 
rippy Power processor efficiency, given by (16) 
rip-j-^ Efficiency given by expression (15) 

ri $ Screen supply efficiency, given by (17) 

Thruster efficiency, given by (11) 
n u Propellant utilization efficiency, given by (13) and (14) 

n u 8r \/ 9l B ’ given (46) and ( 47 ) 

0 In-plane thrust angle 

9* Travel angle increment 

Travel angle 


A 


Primer vector (adjoint to spacecraft velocity) 


A 


X 


c 


x 
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X 


s 


x v 


I 


X 


X 


X 


V 


X 

ir 


A 

T 



P 


y t 


Av 


IT 


IT. 


Magnitude of A ; a general adjoint variable; the iterator 
inhibitor 

Adjoint variable associated with jet exhaust speed 

Adjoint variable associated with reference thrust acceleration 

Adjoint variable associated with degradation time 

Adjoint variable associated with beam voltage Vj 

Thrust cone angle Lagrange multiplier 

Adjoint variable associated with mass ratio 

Adjoint variable associated v/ith performance index (symbol k 
in [1]} 

Adjoint variable associated with propulsion time 

V 

Adjoint variable associated with thrust cone angle 

Gravitational constant of the sun; a general gravitational 
constant 

Gravitational constant of the primary target 
Mass ratio 

Mass ratio increment at an intermediate target 
Performance index; ratio of circle circumference to diameter 
Partial derivative of ir with respect to arbitrary variable x 


xxi 


Auxiliary variable used in equations (79) of [1] 

Thrust switch function, given by (4?) 

Special form of thrust switch function, given by equation (186) 
of [1] 

Portion of total thrust switch function, given by (193) of [1] 
Propul si on-corner-proximi ty to! erance- i nterval 
EPS propulsion time 
Characteristic degradation time 

Transformation matrix for rotating from ecliptic to equatorial 
coordinate system 

Thrust cone angle (between thrust and radius) 

Angle between normal to solar arrays and the spacecraft-sun! ine 

Out-of-plane thrust angle 

Longitude of ascending node of an orbit 

Angular position from the ascending node of an orbit to the 
spacecraft; argument of peri focus of an orbit 


1.0 INTRODUCTION 


This document is the first supplement to the currently existing primary 
HILTOP program document (published in December 1974; see reference [1]) and 
describes the modifications and improvements made to the HILTOP electric 
propulsion trajectory optimization computer program up through February 1978. 

A new, more realistic propulsion system model involving the actual 
ion beam current and voltage has been implemented in the program. The 
power processor efficiency, ion thruster efficiency, and thruster specific 
impulse are modeled as variable functions of the (solar array, nuclear, 
or other) power available to the propulsion system. The number of operating 
thrusters are staged, and the beam voltage is selected from a set of five 
(or less) constant voltages, based upon the application of variational 
calculus. The minimum and maximum number of operating thrusters, the 
minimum throttling ratio, and the maximum input power. to an individual 
thruster are specified as input data. The constant beam voltages may 
be optimized individually or collectively. 

The new propulsion system logic is activated by a single program 
input key (NAMELIST input "NEW"); program modifications have been designed 
to retain the “old" HILTOP program within the framework of the new logic, 
so that old input data files (with no modifications required) will run the 
new program version and produce identical results as before. 

The capability of simulating solar array degradation with the new 
spacecraft model is not included in this program version; also not in- 
cluded is the capability of simulating the new spacecraft model under 
the Launch Vehicle Independent (LVI) mode. The simulation of array 
degradation and the LVI mode remain available with the old spacecraft 
model . 

The execution step requirements of the new program version are a little 
less than 390K bytes of Main Core Storage. This compares to 350K for the 
old version. 
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The report contains the new analysis describing these features, a 
complete description of program input quantities, and sample cases of 
computer output illustrating the new program capabilities. 
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2.0 FORMULATION 


2.1 Spacecraft and Trajectory Models 

The following discussion is oriented toward the programming logic 
aspects of the new HILTOP computer program version. Equations and analysis 
which have not been affected by the implementation of the new spacecraft 
model are not repeated here and may be found in reference [1]. The new 
spacecraft model was obtained from the Lewis Research Center [2] . 

2.1.1 Spacecraft Mass Components 

In the new spacecraft model, the spacecraft is composed of an 
electric propulsion system and associated tankage and propellant masses, 
a structure mass component, a retro propulsion component (for maneuvers 
about a primary target), a set of instrument package masses to be dropped 
at intermediate targets and a net spacecraft mass as follows: 


n-1 

m o = m a + ra th + m p + B t + m s + "V * 2] "drop 1 + m net* 0) 

i=l 

where m Q is the initial spacecraft mass; m a , m p , m^ are the solar 

array (or other power source) , thruster subsystem, propellant and tankage 
masses, respectively; m s is the structure component; m^, is the retro 
propulsion mass; m drop .. is the instrument package mass left at the 
target; and m pet is net spacecraft mass (payload). In the analysis to 
follow, the subscript o denotes the launch body and n the primary (final) 
target. The net spacecraft mass consists of the scientific instruments, 
communi cations, navigation, and other engineering hardware, shielding, and 
any other mass components required to carry out the mission of interest. 
Equation (1) is identical to that of Reference [1] (the old spacecraft 
model) except that the quantity m ps of the old model has been replaced 
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The solar array and thruster subsystem masses are given by. 


m a = Vao + m b 
H'th = a t n max p to 


( 2 ) 

(3) 


where p is the reference power of the solar array or other power source 

dO 

(see Electric Propulsion System), a „ is the specific mass of the arrays, n m _ v 
is the maximum number of operating thrusters, p^. Q is the reference power of 
each thruster, is the specific mass of the thruster and power conditioning 
subsystem, and m^ is a constant mass. 


The propellant, tankage, and structure masses m p , m^, and m s , respect- 
ively, are computed the same as in the old model. The retro propulsion 
mass m r is composed of the two components (as before). 


"V * % + 


m 


rst 


(4) 


in which the retro structure and tankage mass component is computed 
identically as before and the retro propellant requirment, m rp , is now 
given by 


"Vp = Mvn - • J t*tX • (5) 

j a » and are input jettison indicators set equal to one if the solar 
array, thruster subsystem and tankage mass components are to be jettisoned 
prior to the retro maneuver and equal to zero otherwise, and j is an 
indicator (1 or 0) for the presence or absence respectively, of the retro 
stage; and e is the retro burn exponential factor discussed in reference [1]. 
j ft is a new indicator, and j^ replaces j ps of the old model. 

2.1.2 Electric Propulsion System 

The electric propulsion system is comprised of the solar array and the 
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thruster subsystem*. The solar array is characterized by a reference power, 
P ao , defined as the power developed at 1 AU from the sun assuming the 
arrays are oriented normal to the sun. The instantaneous power developed 
by the arrays p , takes into account the effects of temperature, the distance 
from the sun *nd the orientation of the arrays with respect to the sun line. 
The instantaneous power developed is expressed 

" •»P«o < 6 > 

where y is a power function which accounts for temperature, distance and 
orientation effects and is computed identically as in reference [1], having 
the same options. 

The thruster subsystem consists of a power distribution system which 
accepts the power delivered by the array and distributes it to a number of 
power processing units (PPU's) each of which is dedicated to a separate 
thruster. Each PPU processes the power input to the unit to deliver to 
the associated thruster the appropriate voltage and current parameters for 
efficient operation. The power output of the power distribution system, 
p^, is modelled 


Pd = ” P dP a < 7 > 

where is the efficiency of the power distribution system, a specified 
constant. This power is distributed evenly among the current number n t 
of the operating PPU/thruster modules. That is, the power input to each 
PPU is 


„ _ p d~ p conv 

* ' !81 

p conv 1S P ower requirement input to the system DC-DC converter to provide 
power to the digital interface units, the power processor low voltage section, 
the thrust system, and the mission module. P conv > which is loosely denoted 
“housekeeping power 1 ' is given by 

A power flow schematic is displayed in Appendix A. 
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- ( p diu* p gy) n t * p mm + p ts 


p conv 


^conv 


(9) 


in which five program-input constants are 

digital interface unit power requirement, 
power processor low voltage input power requirement, 
mission module power requirement, 
thrust subsystem power requirement, 

DC-DC converter efficiency. 

The numerator of equation (9) represents the power output of the DC-DC 
converter. The power output of each PPU/ thruster module is the beam 
power which is related to the input power as follows 


p diu 

P*v = 
p mm 

p ts = 
n conv 


p b = r 'TH n PPU p t 


( 10 ) 


where rij^ and nppy are the efficiencies of the thruster and power processing 
units, respectively. These efficiencies are dependent upon the operating 
conditions, expressed in terms of the beam current I g and the beam voltage 
Vj , of the thruster. The beam voltage is selected at each instant in time 
from one of up to five input discrete values; the selection is made by the 
program as part of the problem solution as discussed in the section, 
Optimality Conditions. The beam current is throttled as necessary to make 
use of the input power , but is subject to the constraint of a maximum 
operating value, I max , which is specified by input. The thruster efficiency 
is modelled empirically 


n TH = n u n 


PTH 


01 ) 


where ccF^ is a thrust reduction factor due to double ions and beam 
divergence 

aF t = C 1 + z zh 02) 

where C-j and Cp are input constants; n u is the propellant utilization 
efficiency which is a function of the beam current 



ri 


u 


I 


m„ + m + I, 


B 

(1-C 3 I b /C 4 ) 


for Ig < 1 amp 


( 14 ) 


where m. is the thruster neutral propellant loss in equivalent amps* m_ is 
* n 

the neutralizer propellant flow rate in equivalent amps, and m^ , m n , C 3 
and are input constants; and Hpy^ is a function of the beam voltage 

y , 

"pth (l+CjKVj+Vg-aVj) + ej + aVj ' 

where £j is the discharge losses in eV/ion, AVj is the discharge voltage 
in volts, V G is the neutralizer to beam coupling potential in volts, and 
C 5 is the ratio of accelerator current to beam current, ej » aVj, Vg 
and Cg are input constants. The power processor efficiency is of the form 
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(16) 


O+CgHVj+Vg-AVj) + £j + aVj 

n PPU = (l+C 5 )(V I +V G -AV I )/n s + (ej+AVjj/np 


where n D is the discharge supply efficiency, an input constant, and n s is 
the screen supply efficiency which is assumed to be a linear function of 
beam current; i.e.. 


n s C 6 + 

For trajectory computation purposes, it is convenient to express the 
thruster subsystem performance in terms of the actual thrust output f^ by 
the n^ operating thrusters 

ft = ^ newtons ^ ( 18 ) 

-3 

^here k = 2.0391x10 ) and the specific impulse I^p 


T _ 

X SP = kg r 


(sec) 


(19) 


o 

where g r = 9.80665 m/sec . Other performance parameters of traditional 
importance include the jet exhaust speed c , the mass flow rate m , the thrust 
acceleration a , and the beam power p^ . These parameters are evaluated 
as follows: 


' 9 r ! SP 

(20) 

■ - V c 

(21) 

“ V” 

(22) 

■ WtC 

(22) 


where m is the instantaneous mass. 


8 


ORIGINAL PAGE I b 
OF POOR QUALITY 


For identification and documentation purposes, it is useful to identify 
reference values of the thruster specific impulse and power. These 
reference values correspond to the maximum permissible value of beam current; 
I , an P arame t er » an d the largest of the discrete values of beam 

voltage, Vj mQX . The reference value of specific impulse, denoted I Spo , 
is evaluated with Equation (19) using I max and Vj max in the Equations (11) - 
(17) and (19). The unit thruster reference power p^. Q is defined to be the 
power allocable to each thruster that is input to the power distribution 
system assuming each thruster is operating at a beam current of I and a 

Hid a 

beam voltage of Vj max . That is. 


. + * V I , 

t0 n pd n PPU|j „ ' 

■ 1 max ,v Imax 


n D (UC 5 )(v Hriax' l ' V G- aV I ) + (C 6' l ' C 7 1 iiiax )tE l' l ' aV I 

n pd VVVmax* 


Thus, if n max denotes the maximum number of operating thrusters permitted, 
then the maximum power output of the arrays that can be utilized by the 
thruster subsystem and other spacecraft modules is 

p = „ p + !convl SS xl 

a raax max t0 n pd 


in which P conv(max) is P conv evaluated with n t = . 

The maximum power input to an individual thruster is computed. 


p max ’ 10 ‘ 3 *W[ (1+C 5 >< v t ♦ <*l +aV I>] 


9 


(in kw) when I max is specified as a program input, and this formula is 
inverted to compute I when p max is specified as a program input. The 
minimum throttling ratio is defined in terms of the beam current, 

Vatio " Vin^max 

when I - n is specified as a program input, and this formula is inverted to 
compute I . when t t .j 0 is specified as a. program input. 

Notice that the solar array tilt angle does not appear in this formulation; 
it is essentially replaced by Ig. When I g = I mQX , which corresponds to 
"operating below the power curve", then this model is not concerned with how 
any potentially-available excess power is avoided. Methods of avoidance 
include 


(1) tilting the (solar) arrays 

(2) shielding the arrays 

(3) dumping excess power via radiators 

(4) shunting the excess power for other use 

The solar arrai tilt angle would become a factor in the model if array 
degradation were considered; however, including array degradation in the 
model would introduce considerable complexity to the equations and 
algorithms required for the optimal solution and is deemed beyond the scope 
of the current implementation. 

When the spacecraft is "operating on the power curve" as specified by 
Y (i.e., using all available input solar power), the beam current for each 
operating thruster is given by 


I 


B 


-b -^/h 2 + 4CgCyGp t 
2C ? G 


(26) 
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in which , using equations (6), (7) and (8), is given by 


n pd p ao Y p conv ■ 

' ” n t n t 

(27) 

G " Uj + AVj)/n D 

(28) 

fc> = Fy ■ C 7 p t + C 6 G 

(29) 

F V = O+CgXVj+Vg-AVj) 

(30) 


In this spacecraft model, the solar array (or other power source) may not 

be perfectly matched to the thrust subsystem capability under reference 

conditions; instead, the power, p , may be a specified constant, which 

may be less than, equal to, or greater than the thrust subsystem maximum 

power requirement p . In the case of solar electric propulsion, the 

max 

analyst pre-selects the desired power curve y via program inputs a. 

(power curve coefficients) and y - In fact, the same power curve options 

are available as in the simpler spacecraft model, and these options are 

discussed on pp. 11-13 of [1). Then, at any given point in time, the value 

of y will allow the computation of a hypothetical value for the beam 

current using equation (26), and if this value lies in the acceptable 

range I. ,< I D $ , it is used in the actual calculations; if the 

value is greater than I , then I max is used, and if less than 1^ 

(which is associated with the minimum throttling ratio), then special 

action is required involving an iteration to isolate the point along the 

trajectory at which I D = I . , at v;hich point the control state of the 

o mi n 

thruster subsystem is optimally switched (as discussed in the section. 
Optimality Conditions). When either I D = I „ or the spacecraft is 

O IHaX 
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operating on a level portion of the power curve v = Y max > this model is 
not concerned with how any potentially-available excess power is avoided 
(as mentioned above). In most applications, the spacecraft will be operat- 
ing on the non-constant portion of the power curve, representing the optimal 
situation in which the arrays are normal to the sun line and gathering as 
much power as possible. If the program user wishes to have the solar array 
output power matched to the remainder of the spacecraft (under the reference 
conditions of having the solar arrays at 1 AU from the sun and oriented 
normal to the sun line), a program input key is provided which causes the 

setting p = p internally, 
max 

2.1.3 Differential Equations 

The differential equations of motion applicable with the new space- 
craft model are (Consult Nomenclature for definitions of previously- 
defined symbols): 


■ ae t ' "7 R • 

(31) 

r 


* _ 

v r * 

(32) 

I = 1,2, 3,. ..(max 5) , 

(33) 

$ = 0 , 

(34) 


in which the thrust acceleration is given by 


a 


KCi+Vb)*! n t 1 B 


m v 
acc o 


( 35 ) 


and the jet exhaust speed is given by 
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c 


2( c 1 +C 2 I B )n l /^ 


(36) 


where c acc = 5.9301282604x10" ^ m/sec^ and c- ve -j = 29784.916613 m/sec 
render the quantities a and c expressed in program internal units. It is 
therefore emphasized that the symbols a and c from this point of the 
discussion onward pertain to thrust acceleration and jet exhaust speed 
expressed in (normalized) program internal units, in contrast to the a 
and c of equations (22) and (20), respectively, which are expressed in 
MKS units. Relations (35) and (36) are valid regardless of the algorithm 
by which the beam current Ig is generated. The controls are n.., Vp , 
and 1^ , where n^ is the number of operating thrusters, Vj is the beam 
voltage, e^ is a unit vector defining the direction of thrust, and Ig is 
the beam current per thruster. 

The analysis pertaining to the thrust cone angle $ (the angle between 
the radius and thrust vectors) remains essentially unchanged compared to 
that of the old HILTOP model, and is included here for the sake of 
completeness. The propulsion time t of the old model is not included in 
the new model since a spacecraft having thrusters which are staged accord- 
ing to available power does not have a "propulsion time" which can be 
simply implemented in a variational calculus approach. 

The differential equations which govern the behavior of the adjoint 
variables are given by 


.. o * ah T n 

A = — Tjr A + — g- (A*R)R + -gg- + - r eos$) , 


(37) 


X = - (A--S1) , 

V V t' 


(38) 
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(A*e t ) 

Vl^r + 


(A-e t ) 


C 1 +C 2 I B, 


I 

+ ((A'? t )-a) ~ +a 


31^ 

av: 


I = 1 ,2,3,. . . (max 5) , 


( 39 ) 


^ = \ R * ( mxe t) (40) 

The term hj in equation (37) is the component of the variational hamiltonian 
containing the thrust acceleration: 

h r = aa , (41) 

where the coefficient a of the thrust acceleration is called the thrust 
switch function and is given by 

cf = A-e t - ^ (42) 

The partial derivative ahj/aR is a somewhat lengthy expression, determined 
as follows; it is first written as the product of partial derivatives 


3hj 3hj 3lg 

3R~ = 31^ "iR 


(43) 


When Ig = 0 (during coast), = 0 and also ahj/3R = 0 because 3lg/3R = 0 . 
Therefore Ig > 0 in what follows. The following partial derivatives are 

ORIGINAL 
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determined in a straightforward manner and the results are given: 





C 1 +C 2 I B 



(44) 


in which the "singularity" 


1 

a 


is removed by using 


hj/a 


a 


Then 


3lg " ( A * e t“ 0 ) 


C 2 , ^ 

C f C 2 l B \ 


(45) 


and 



when I B > 1 ; g = (I B -1) 2 and g* = f&- = 2(I B -1) . 

B 

The partial derivative shj/aig is then determined. It remains to determine 
3l B /aR . When I g = 0 or I g s I max , 3l g /aR a 0 . It therefore remains to 

determine aI g /aR when the spacecraft is operating on the power-constraint 
curve (using all available input solar power). For this case, I g is computed 
from equation (26), and, after some manipulation, 
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( 48 ) 


2C,G-b 

1 + — 2 

2C ? GI B +b 


in which 


3 °t _ n pd^ao dy 
3R n t 3R 


and 


_ .. _2_ -r lx 
3R r 3 e r 3 d 


d l£-± !?t 

dR 2G SR 


(49) 


(50) 


In equation (50), is the radius unit vector from the sun to the spacecraft 
and 3y/ad is generated the same as in the old model. 


In equation (39), 3lg/3Vj 


is given by 


31 

3V 



- <1*C 5 } 




2CyGlg*rb 


(51) 


when the spacecraft is operating "on the power curve", and 3Ig/3Vj = 0 when 
*B = ^ax or X B“ 0 * 


The expression for X in equation (37) is identical to that of the old 

A 

model when written in terms of the thrust acceleration. 


A* (nixe*) 

X = - a — (52) 
R*(mxe t ) 
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if the thrust cone angle is held constant and X is identically zero if 

A 

the cone angle is unconstrained. The unit vector in' is defined the same 
as in the old model. 


m = 


Rxa 
[ Rxa | 


( 53 ) 


From equation (39), it is seen that five (or less) additional differen- 
tial equations (for Xy^) are integrated along each trajectory, but only 

when at least one beam voltage transversal ity condition requires satisfaction, 

as requested by user program input. Actually, each x w can be non-zero 

V I 

only when the corresponding beam voltage Vj is being used to drive the space- 
craft, and is identically-zero otherwise; therefore, at any given time 
along a trajectory, only one Xy need be integrated. (Transversal ity is 
discussed in a later section.) 


2.1.4 Optimality Conditions 

The control variables available for optimization along the trajectory 
consist of the number of operating thrusters n^, the beam voltage Vj 
(selectable from a set of from one to five constant values), the unit vector 
e^ defining the direction of thrust, and the beam current per thruster Ig . 
The control variable h a of the old model is replaced by Ig , such that 
Ig h 0 defines a coast phase. As will be explained below, Ig = 0 is used 
to conceptually define when the spacecraft is coasting rather than 

= 0 or Vj = 0 , since the optimal values of n^ * 0 and Vj f 0 must be 
maintained even during coast. (Of course, the voltage can be turned off 
on board the spacecraft during coast.) 

The application of the Maximum Principle of optimal control theory 
leads to the result that the proper choice of the control variables is 
that which maximizes the variational hamiltonian, h y , at each point 
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along the path. The variational hamiltonian for the problem forum! ated 
here may be written 


h y = ao - (A-R) - A-R (54) 

r 

where a is the thrust acceleration given by equation (35) and a is the 
thrust switch function given by (42); a gets its name "thrust switch 
function 1 ' because its sign determines, as usual, if the spacecraft is 
thrusting or coasting. Specifically, the variational hamiltonian is 
maximized with respect to the beam current Ig by setting Ig to its lowest 
permissible value (zero) when a<0 and its largest permissible value when 
o>0 , for given values of and Vj . The method of selecting and main- 
taining the optimal values of n^ and Vj is discussed below. For now, 
assume that the optimal values of n t and Vj are known. Then, when a>0, 

Ig is chosen as large as possible. This is accomplished by first 
computing the beam current according to (26), and using this value for Ig 
unless it is greater than I max » in which case Ig = I max is chosen. When 
cr<0 , Ig = 0 is the optimal choice. In this manner, the optimal choice 

of I R is determined. There is also a value of the beam current, I . , 
d mm 

determined by the minimum throttling ratio, which does not enter into the 
discussion here for choosing the optimal beam current, but which plays a 
fundamental role in the algorithm for choosing optimal values for n^. and 
Vj . This will enter into the discussion below. 

The variational hamiltonian h y is maximized with respect to e^. 
identically the same whether employing the old or new spacecraft model; 
the optimal choice for is discussed on p. 16 of [1]. 

It remains to determine the choices for n^ and Vj which maximize h y 
at each point along the path. In general, n^ and Vj might be considered 
to generate a 2-dimensional control sub-space in which is found a finite, 
bounded, discrete rectangular set of permissible points: 
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V 


I 


max 






• • 



mm 


n„. 

min 


■ 9 


9 

-> 


9 9 

^ax 


Then, at the initial point in time, the variational hamiltonian is computed 

for each point of the entire grid (for each pair of values (V^n^)), and 

such that Ig is chosen as a function of Vj and n^ according to the method 

discussed above* Cases for which I D < I . are discarded from the 

fcs min 

competition. When the grid mapping is complete, the optimal values of n^. 
and Vj are known, by simply saving the values associated with the largest 
h y as the grid mapping progresses. It is assumed in the present 
implementation that all cases are not discarded from the competition for 
maximum h y due to Ig < I . n# i.e., there is sufficient power available to 
thrust at the start of the mission for a sensibly designed spacecraft. The 
optimal values of n^ and Vj thus determined might be represented by the 
circled point in the grid above. If the thrust switch function ct associated 
with these optimal values is positive, the spacecraft thrusts (using the 
maximum value of Ig allowed); otherwise, a<0 and the spacecraft coasts. 

* ■ 

This is depicted in Appendix B. 
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A critical distinction regarding the beam current Ig is made at" 
this point. The notion of a "phantom" beam current Igp is introduced 
and is defined as the value of beam current (per thruster) which the space- 
craft would have if there were no physical restrictions (I . n *I raax ) on 
the beam current. Igp is therefore computed using expression (26). Also 
introduced is the restricted phantom beam current, Ig X , which is equal to 
Igp unless Igp > I max , in which case Ig X = I max . Therefore, both Ig p and 

I DV may be less than I . , but only I Dn may be greater than I. The 

da min Dr JTiaX 

actual beam current, which is used to drive the spacecraft when a>0 , 

cannot violate the (throttling ratio) bounds I « n ? Ig $ I max except for 

I B = 0 . Both r BP and Ig^ are always greater than zero. 

It is the restricted phantom beam current Igj, which is always employed 
in the computation of the jet exhaust speed, c, in equation (36), even 
when Ig = 0 . The jet exhaust speed (or equivalently, specific impulse) is 
therefore always a positive quantity, even when the spacecraft is not 
thrusting. The specific impulse is therefore conceived as a latent 
physical property which characterizes a thrust subsystem's potential 
capability, even when the spacecraft is not operating. In turn, t* n jet 
exhaust speed c, always computed using I gx , is used in equation (42) for 
the thrust switch function a. Then the sign of a is used, in turn, to 
decide whether the spacecraft is thrusting or coasting, in which a is 
computed using the optimal values of Vj and n t . Therefore, using functional 
notation. 


I 


BX 



(55) 


Then, in terms of Vj and n t , the jet exhaust speed is given by (again 
using functional notation) 


C = (c V g 1 k)( C l +C 2 I BX^ V I ,n t^) n u( I BX^ V rV) ,/ ^7 


(56) 
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Once optimal values for Vj and n t are initialized, it remains to 
determine how to maintain optimal values along the path. A brute force 
approach would be to keep testing the entire grid of permissible points 
in the (V^n^) control subspace along the path, and switch to a new 
value of (Vj,n t ) when the new controls produced a larger variational 
hamiltonian; the point of switchover would be isolated by iteration. 
However, for control grids having 30 or 40 points, this approach would be 
computationally very costly. An assumption is therefore made which is 
denoted "the neighboring solution assumption," for the purpose of building 
a computationally efficient optimal control algorithm. In this assumption, 
the only grid points in the control subspace which are considered as 
candidates for optimal control are those "directly" neighboring the 
presently existing optimal control point; the neighboring points are 
effectively defined by the connecting lines (to the circled optimal point) 
in the grid depicted above. In a two-dimensional control grid there are 
therefore (a maximum of) eight alternate control strategies to consider, 
and, obviously, fewer when the optimal control point is on a grid border 
or corner. 

The tactic employed in the computer program for maintaining optimal 
values of Vj and n t therefore consists of the above scheme of comparing 
candidate hamiltonian values, in which each point (in time) of optimal 
switchover of the controls V- and n^. (i.e., the point where the difference 
between hamiltonian values associated with the optimal grid point and 
candidate grid point vanishes) is strongly isolated by iteration. 

A few other considerations must be taken into account in the algorithm 
for maintaining optimal values of Vj and n^ . Roots of the two functions 
I Bp ** I max -and Ig p - 1^ must be strongly isolated (by iteration), 

corresponding respectively to the maximum throttling ratio (unity) and 
minimum throttling ratio thresholds. It is useful to introduce at this 
point the notion of an "imaginary spacecraft" associated with each point 
in the (Vj, n^) control subspace neighboring the optimal point, such that 
the imaginary spacecraft has the corresponding neighboring values of beam 
voltage Vj and number of operating thrusters n^. . Then minimum throttling 
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ratio thresholds must also be isolated for the imaginary spacecrafts 
associated with {some of the) points neighboring the optimal point in the 
(Vj, n^) control subspace, because switches in this subspace are not allowed 
which would result in I g < I .j n . Switches in the (Vj, n^) control subspace 
will occur regardless of the sign of the thrust switch function a . 

When the minimum throttling ratio threshold is attained and I^p is 
negative, a switch must occur in the (Vj, n^) subspace or else the space- 
craft must commence coasting unless it was already coasting (a<0). Or, if 
the throttling ratio of any comparative imaginary spacecraft attains the 
minimum threshold and Ig p for that imaginary spacecraft is positive, then 
a switch in the (Vj, n^) control subspace may occur to the point associated 
with that imaginary spacecraft. It may also be that, for solar electric 
propulsion, the spacecraft will recede so far from the sun that all points 
of the (Vj, n^) control subspace grid will "fae u below the minimum 
throttling ratio threshold, so that the spacecraft will have no alternative 
but to coast; and it may also be that the spacecraft will commence thrusting 
again as it approaches the sun and some pair of values (Vj, n^) have an 
associated acceptable throttling ratio with a>0. Whenever one of the 
above situations (described in this paragraph) occurs, the primer 
derivative must be jumped so that the variational hamiltonian remains 
constant. The jump condition is 


* i 

= A 



R 


(57) 


in which hj = a a is that portion of the variational hamiltonian associated 
with the engine parameters. 

It is possible that the neighboring solution assumption (discussed 
above) will in fact be violated at some point in time along a particular 
trajectory, for a particular spacecraft configuration and a particular 

mission. This situation will occur extremely infrequently. However, when 
it does happen, the optimal controls (Vj,n t ) must be re-determined. It is 

- gSES® 


a simple matter to detect when the neighboring solution assumption has been 
violated; whenever a switch occurs in the (Vj,n t ) control subspace to a 
point neighboring the prior point, the values of all variational hamiltonians 
associated with the points neighboring the new optimal (?) point are compared 
to the hamiltonian value associated with the new point. If the new point has 
the largest hamiltonian, it is (considered to be) the optimal point (Vj,n^.), 
and the neighboring solution assumption is not violated; otherwise, the 
assumption is violated, and one of the new neighboring points in the (Vj,n t ) 
control subspace is more optimal (the one having the largest hamiltonian), 
and an additional switch is immediately made to the new point, with a jump 
in A according to (57) to maintain hamiltonian constancy. Once the switch 
to the new point has occurred, the entire test of neighboring points for the 
optimality of the present values of (Vj,n^) is conducted once again. The test 
of the neighboring (optimal) solution assumption is indeed carried out every 
time a switch occurs in the (Vj,n^) control subspace, and the testing of 
neighboring points is repeated (as described above) until the optimal values 
of (Vj,n^) are obtained. 

Whenever the neighboring solution assumption is violated (which, again, 
happens extremely infrequently), a double-switch or multiple-switch occurs 
in the (Vj,n t ) control subspace, with a concurrent jump in A ; this represents 
a slightly sub-optimal control strategy, compared to the (computationally 
inefficient) globally optimum strategy of directly testing for the optimality 
of all points in the (Vj,n t ) subspace. The globally optimum strategy 
would therefore switch directly to the (not-necessarily-neighboring) new 
optimal point at a time along the trajectory very slightly less than that 
found by the neighboring solution strategy, and with no jump in A required. 
Nevertheless, the mechanization adopted for the general solution to the 
optimal rocket' flight problem, as posed in this report, employs the neighbor- 
ing solution assumption (or restriction) because: 

» It is computationally efficient 

* It is extremely rarely violated 

a All violations are felt to have negligible impact on spacecraft 
masses and other performance parameters (compared to the globally 
optimum solution) 
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2.2 Boundary and Transversal ity Conditions 


The basic boundary conditions and transversal ity conditions are described 
in [1], This section discusses only those conditions which must be modified 
to accomodate the more sophisticated spacecraft model, and also some entirely 
new transversal ity conditions. Therefore, any boundary or transversal ity 
conditions appearing in [1] but not mentioned here remain unchanged. 

The discussion in this section pertains solely to the new spacecraft 
model, and therefore certain quantities which are indigenous to the old 
model will be found entirely absent in the new analysis. Specifically, 
these quantities are the (constant) jet exhaust speed c of the old model, 
the reference thrust acceleration g, and the propulsion time t. Therefore, 
all equations in [1] pertaining to these quantities are absent in the 
analysis describing the new spacecraft model, and any boundary or trans- 
versal ity conditions associated with these quantities are not applicable 
with the new model . 

The general equation for the transversal ity conditions is written 
X dir + 

IT 


Ilitljiax 


■2 [ A * dR ~ A ’ dR + \> dv + y* % w h v dt -r. 

i=i ■ i=i 1 i-i 


It] 1 ' -0 


(58) 


in which .he performance index ir is equal to the negative of the net space- 
craft mass, t r = - m net , and is the arbitrary positive constant which 
renders the general transversal ity condition linear and homogeneous in the 
adjoint variables (x^ replaces the symbol k in [1]). The convenient choice 
is made whereby each is forced to be continuous at each intermediate 
target, which means that only Xyj(t n ) need a PP ear in derived trans- 
versal ity expressions rather than the cumbersome expression 


(59) 





I - 1 ,2,3,.. .(max 5) 


This is because Ay^tj,) alone, with A Vl (t Q ) = D and Ay^{t.) = Ay (t..) for 

each i, has the same value as the cumbersome expression cited above if 
Ay (t Q ) were not zero and Ayjft.j ) were not continuous, and this is due to the 
absence of Ay^ in the differential equations. 


The expression for ir may be written 


* * j r"Vs + ra of k s +k t ' (1+k t>V J V (1+k rt )e x[ (1+j t k t>V Vt 0 + 

n-1 _ n-1 n-1 

1 | 1 ( k samp r k drop d)j + 0 + k t ) ^ampf " k t ^ k drop d 


+ ^ao' l Vf 1 - J aV 1+k rt >e x] + KWto 1 r (1+k rt )e xJ ■ (60) 


where symbol definitions may be found in Nomenclature, it may be written 
functionally in its most general form. 


* = ’ tv - 0 - v «n ,v n* \, x >6 - 1( 


(61) 


Using the notation n x = 8ir/ax , the general variation of it may be written 


dir ® m , dv + ir.. dv + it dv„ + it,, dV T -Hudfi + n.di (62) 

Vcb 0 “O v„ n -n v n n Vj Imax & i 


*max 
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Then, the following consists of a summary of those transversal ity 
conditions published in [1] which are altered when employing the new 
spacecraft model . 

Equation (98) of (1] for Tr m remains the same except the two terms 
containing j are eliminated. Equation (100) of [1] for ir v ren , a -j ns 

unchanged except the term -jp S m ps is replaced by -J a m a - j^th * Further- 
more, there is a transcription error in that equation, such that the symbol 

v should be replaced by v . 
c p 

As before, the arbitrary positive quantity X^ is assigned a value 
which causes the transversal ity condition associated with the final mass 
ratio to be satisfied: 


K = “ K ^ 

n n 

This is defined as the negative of the symbol k in equation (144) of [1]; 
this sign is merely a matter of convention, so that the quantity (-k) 
is simply replaced by X„ in the transversal ity equations appearing in [1]. 

(The form of x^ expressed by (63) above is coded in the program.) 

For open launch excess speed with m Q being independent of <5 and i, 
the transversal ity condition is 


X it 
Ti V 




1 = 0 . 


(64) 


For cases in which m is a function of 5 and/or i, the transversal ity 

o 

conditions are written as follows: for open excess speed. 


“mo 

sv .o 



(65) 
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For open geocentric declination of V 


»q 


ani r 

X IT -r-r 2 - - A * (V xriJxV /V COS6 I - 0 

IT m Q 86 0 L' «Q p' «o wo j 


( 66 ) 


In addition, .the .factor f is dropped from equation (151) of [1], No 
changes were required in the software for the conditions discussed in this 
paragraph, since they are coded in the form given here. 


New transversal ity conditions arise due to the presence of the beam 
voltages Vj in the general transversal ity condition (58). For voltages Vj 

the transversal ity conditions (whose satis- 


less than the maximum V 


Imax 


faction generate optimal values of V j ) are, simply. 




I = 1,2, 3 : 


(67) 


The constant beam voltages may be optimized individually or collectively, 
as specified by user program input. They may also be driven to specific 
values by the iterator. 

A cautionary note is issued at this point to the program user wishing 
to optimize any or all of tne beam voltages. Specifically, the quantities 
X\f j. remain quite small (relative to unity) over any trajectory, 
so that their final values (equation (67)) are inherently small. This is 
due to the fact that =0, I = 1,2,3,..., and the derivatives X 

are of the order of Vj 1 in magnitude. Furthermore, the independent variables 
Vj are relatively large compared to other independent variables such as A q 
and A q , which are of the order unity. Consequently, when attempting to 
optimize the beam voltages, the user must set the corresponding iterator 
independent variable weights to (e.g., ) Xi(5) = 10“^ and dependent variable 
tolerances to (e.g.,) Yi(3) « 10 . However, since only limited experience 

has been gained in optimizing the beam voltages to date, the choice of best 
values for these weights and tolerances is not well understood. 
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Since V Imax appears in the expression for the performance index * 
(through p t ), its transversal ity condition is given by 


X /V T + * v (t ) = 0 
max ^ ax 


(68) 


When 


the power source is not matched to the thrust subsystem, 


9 p^. 

*V T = a t n max^ ~%V 1+k rt) e xW 

W max 


(69) 


When the power source 

where p is given by (25)) 
a max 


is matched to the thrust subsystem (i.e. s P a0 = P^. 


"V = \ax{ ( V^>A (1+ VtK ( Va +j th“t>} W 


8 Pto 


(70) 


max 


max 


In (69) and (70), 


9 Pto _ 


3V, 


X max 


(i+c 5 )i max 

V^ C 6 + ^7 X ma?" 


P t0 C ? (1 +C 5 ) ( V l ma x* V R ~ AV I ^ X max 


max 


T1 pd^ C 6 +C 7 I max^ 


al 


max 


3V, 


(71) 


max 


When I max is specified as constant. 


al 

TV 


max_ 5 0 , 


(72) 


max 


and when the maximum power Input to an individual thruster is specified 
via (25a), 
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( 73 ) 


81 

8V 


max 

^ax 


2 


( HC 5 }I max 


2.3 Auxiliary Computations 

This section presents equations employed in computations which are made 
after the iteration sequence involving the primary target is completed. 

2.3.1 Additional Block Print Variables . 

A standard print block is employed for printing information at various 
points along a trajectory. Each standard block contains a total of forty 
parameters , which are described in [!]• 

The standard block may now be augmented in a third way, in addition to 
the two ways (power degradation and target-relative coordinates) described 
in [1]- When the improved spacecraft model is invoked by program input 
NEW, two additional lines automatically appear as the sixth and seventh lines 
of the print block. The information contained in these lines is as follows: 

NO. THR Number of operating thrusters, n t . 

VOLTAGE Beam voltage, Vj, in volts. 

CURRENT Beam current per thruster, Ig , in amps. 

PHAN CUR Beam current per thruster I Bp which would 

be realized if there were no limits imposed on 
the beam current; in amps. ("Phantom current", 
given by expression (26)). 

SP IMP Specific impulse, I $ p, given by expression (19), 

in seconds. 
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THRUST 
BEAM POWER 

DUMP POWER 

UTIL EF 

THR RED 

THR EF 

PPU EF 

PTH EF 
SS EF 

ARRAY POWER 
PPU POWER 


Total thrust, f^ , as given by expression (18), 
but expressed in pounds. 

Beam power per thruster, as given by expression (23) 
except divided by the number of operating thrusters 
n t , in kilowatts. 

Dumped power, in kilowatts; zero when I D $ I , 

D FTlctX 

otherwise computed from v j ^ I BP“ I max^ * where j bp 
is the phantom current defined above. 

Propellant utilization efficiency, n u , as given 
by expression (13) or (14). 

Thrust reduction factor, aF t , as given by 
expression (12). 

Thruster efficiency, n TH , as given by expression 

( 11 ). 

Power processor efficiency, Tippy , as given by 
expression (16). 

Efficiency np^ given by expression (15). 

Screen supply efficiency, n s , given by expression (17). 

Total po wer output by the solar arrays (or other 
power source), p a , as given by expression (6), in 
kilowatts. 

Power input to each PPU module, p t , as given by 
expression (8), in kilowatts. 
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Whenever one of the triggers associated with the iterator independent 
variables X36 through X40 (the beam voltages) is turned on* the 
corresponding adjoint variables are integrated and printed as an 
additional line of the print block following the output quantities 
described above. 

2.3.2 Additional Extremum Point Summary Print 

When the improved spacecraft model is invoked by program input NEW, a 
second table of values is printed beneath the currently existing "Extremum 
Points of Selected Functions." 

The usual extremum table contains a cross-reference column at the far 
right, which contains summary information relating to the second table if 
the entry (line) was caused by a parameter in the second table; otherwise 
the cross reference field is blank. The cross reference information con- 
sists of IMAX to denote points at which I n attains the I mav threshold, IMIN 
to denote points at which Ig attains the 1^ threshold, +_ V to denote beam 
voltage switch points, and +_ N to denote points of thruster staging. 

The second table contains the time (repeated from the first table); 
the number of operating thrusters; the beam voltage (volts); the beam current 
per thruster Ig and phantom current Igp (defined in the preceding section), 
both in amps; the propellant utilization efficiency n u ; the specific 
impulse I^p in seconds; the beam power per thruster p^ in kilowatts; the 
dumped power in kilowatts, as defined in the section immediately preceding, 
and the total thrust f^ in pounds. 

If the engine state (number of thrusters and/or beam voltage) has 
switched, the corresponding table entry will be flagged with a plus sign or 
a minus sign, and the entire line will be repeated with values corresponding 
to after the switch has taken place. The plus and minus signs will also 
appear during coast phases, when the number of thrusters and beam voltage 
are printed as zero. Also, critical values of beam current are flagged 
with an asterisk(*). 


3.0 PROGRAM INPUT 


The following consists of a complete description of program inputs. 
With respect to the basic HILTOP report [1], many new input quantities 
necessary to characterize the more realistic propulsion system have been 
added. All new input names are flagged with a single asterisk (*), and 
inputs whose definitions have been modified are flagged with a double 
asterisk(**). 

The new propulsion system logic is activated by a single program input 
key (NAMELIST input "NEW"); program modifications have been designed to 
retain the "old" HILTOP program within the framework of the new logic, so 
that old input data files (with no modifications required) will run the 
new program version and produce identical results as before. 

3.1 Namelist 


Inputs to HILTOP are given through the NAMELIST feature of the Fortran 
programming language. The input NAMELIST is named MINPUT, and every input 
required or used in the program is declared by name in the list. The general 
form for assigning an input value to a quantity is, simply, 

NAME=VALUE 

where NAME is the name assigned to the variable and is included in the NAME- 
LIST, and VALUE Is a numerical or logical quantity consistent in form (i.e., 
logical, integer, or real) with NAME. Unless otherwise specified, all 
MINPUT names commencing with one of the letters T through N represent 
integers, whereas all names commencing with one of the letters A through 
H or 0 through Z are double precision floating point numbers. Each 
NAMELIST case must begin with the characters 

BLANK. f'&JT* 

&MINPUT 

commencing in card column 2 and followed by at least one blank, and end with 
the characters 


33 


&END 


preceded by at least one blank. Card column 1 is ignored on all NAMELIST 
input cards. Multiple data assignments on a single card are permissible if 
separated by commas. Blanks in the variable field, VALUE, are taken as 
zeroes. A comma following the last VALUE on a card is optional on the IBM 
system. The order of the input data assignments is arbitrary; i.e., they 
need not be in the same order as listed in the NAMELIST. In fact, there is 
no requirement that any specific input parameter be represented in the input 
data set. If no value is included in the inputs for a particular parameter, 
the default value is used (see Default Values). For other details regarding 
the NAMELIST feature, the reader is referred to, for example, the IBM 
System 360/Fortran IV Language manual. NAMELIST cases may be stacked 
back-to-back indefinitely. A single NAMELIST input error may cause the 
remaining NAMELIST inputs to be ignored. 

3.2 Definitions of Input Parameters 

Specific examples of the program inputs are given in the Sample 
Problems and Results section. Default- values of inputs are given in the 
next section. 


The program inputs, in alphabetical order, are: 


AAI Desired final extra-ecliptic inclination, i. Related to AE, 

AR, and IOUT. [deg] 

AE Desired final extra-ecliptic eccentricity, e. Related to AAI, 

AR, and IOUT. 

**ALPHAA Specific mass of solar arrays, a a . (See expression (2) of 
this report or expression (4) of [1].) [kg/kw] 


**ALPHAT Specific mass of power distribution, processing, and thruster 

subsystem, a t . (See expression (3) of this report or expression 
(4) of [!].) [kg/kw ] 


ALTITU This input variable is associated with program logic which has 

not been kept up-to-date, specifically, logic pertaining to opti- 
mum departure of a NERVA-type rocket from earth orbit. This 
variable should be ignored. 
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AN 

AR 

ASOL 


BI 

*BMASS 



CNI 

CNIX 


*CSEP 

CSTR 

C'TANK 

CTRET 

*CURMAX 

*CURMIN 

*CVOLT 


Trajectory- integration exponent n in expression (39) of [1]. 

Desired final extra-ecliptic perihelion distance, r^ . Related 
to AAI> AE, and IOUT. [AU] 

Array of five elements consisting of the solar power law co- 
efficients a- in expression (18) of [1] • AS0L(1)^0 tells 

the program to use the input coefficients rather than the internal 
coefficients. The coefficients are normalized internally, and 
the program executes the iterations to produce the required 
remarkable points of the power curve (which are printed). 

Efficiency coefficient b in expression (16) of [1]. Related 
to DI and El. 

Constant nass in expression (2), m^ . [kg] 

Launch vehicle coefficients b-j, bg, and b 3 in expression (2) of [1). 
Used only if MBOOST is negative. 

[kg, m/sec, kg] 

Inclination to ecliptic of primary- target orbit. Input only when 
M0PT3 = 11. Related to ECI, OMI, SAI, SOI, TPI, EMUODD, and 
RADODD. [deg] 

Array of five elements, the first three of which may be currently 
used. Inclinations to ecliptic of intermediate-target orbits. 

Input CNIX(i) only when MOPTX(i) = 11. Related to ECIX, OMIX, 

SAIX, SOIX, TPIX, EMUODX, and RADODX. [deg] 

Array of seven elements, consisting of the quantities C-j, Cg, C 3 , 
C 4» C 5 S C 6* and *7 9 res P ect i ve ty> found in equations (12) through 
(17). 

Structural factor, k s , in expression (8) of [7]. 

Propellant tankage factor, k^., in expression (7) of [1]. 

Retro tankage factor, k rt , in expression (11) of [1]. 

Maximum allowable beam current for an individual thruster, I 
Related to POWMAX. [amps] max 

Minimum allowable beam current for an individual thruster, I . . 
Related to TRATIO. [amps] min 

Neutralizer to beam coupling potential, V r , in expressions (15) 
and (16) [volts] 
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*DEFFIC 

DI 

*DLOSS 

DMRETR 

**DPOW 


*DVOLT 

ECI 

ECIX 

El 

EMUODD 

EMUODX 

*ETAPD 

*ETCONV 

GAMMAX 


Discharge supply efficiency, np » in expression (16). 

Efficiency coefficient d in expression (16) of fl]- 
Related to BI and El. [km/sec] 

Discharge losses, ej , in expressions (15) and (16). [eV/ion] 
Retro engine mass, m rS , in expression (11) of [1]. [kg) 

Ratio of housekeeping power p^ to reference power p .p Used 

only with the old spacecraft model. The power transmitted to 
the propulsion system is that generated by the arrays less 
housekeeping power which is constant along the trajectory. 

The power output of the arrays normal to the sun at 1 AU is 
p re f+ P h . This option should not be invoked on missions 

during which large solar distances are encountered where the 
power developed is less than p^ . Erroneous results will be 
obtained. 

Discharge voltage, AVj, in expressions (15) and (16). [volts] 

Eccentricity of primary-target orbit. Must be less than unity. 
Input only when M0PT3 =11. Related to CMI, OMI, SAI, SOI, 

TPI, EMMUODD, and RADODD. 

Array of five elements, the first three of which may be currently 
used. Eccentricities of intermediate-target orbits. Input 
ECIX ( i ) only when MOPTX(i) = 11. Related to CNIX, OMIX, SAIX, 
SOIX, TPIX, EMUODX, and RADODX. 

Efficiency coefficient e in expression (16) of [1]. Related 
to BI and Dl. 

Gravitational constant of primary-target. Input only when 
M0PT3 =11. Related to ECI, CNX, OMI, SAI, SOI, TPI, and 
RADODD. [m 3 /sec 2 ] 

Array of five elements pertaining to the gravitational constants 
of intermediate- targets. These inputs must be ignored at 
present. 

Efficiency of the power distribution system. n oc p in 
expression (7). . 

DC-DC converter efficiency, n conv > in expression (9). 

Maximum permissible value of the power function y when mode = 5. 
At solar distances less than the value for which y = GAMMAX, the 
solar arrays are assumed to be tilted or shielded such tnat y 
is maintained at the limiting value. 
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**GAP 


HOUR 

IBAL 

INTPR 

IOUT 


IRK 

IRL 


Propulsion-corner proximity tolerance- interval , a a. See 
discussion in the section Avoiding Corners *n the Propulsion- 
time Function in [1]. Whenever the thrust switch function a 
grazes the zero-axis within the tolerance 1 4cr | on any 
trajectory, an internal counter is incremented, and the 
trajectory is considered to be in the neighborhood of a 
propulsion-time corner. Positive value of GAP causes 
forced-thrusting case to be inserted, negative value causes 
bypass to next case, whenever the internal counter reaches 
the related input variable NHUNG. Value is set negative when 
new spacecraft model is invoked. 

Hour-of-day of reference date (e.g., 17.352D0). Related to 
MYEAR, MONTH, and MDAY. 

Ballistic option indicator. Setting IBAL r 0 invokes option 1 
discussed in the section Ballistic Trajectory Option of [!]• 

Indicator which specifies print-length when the iteration in 
subroutine INTERP, fails. Value of 0 causes shortprint and 1 
causes detailed-print. 

Extra-ecliptic mission indicator. IOUT = 1 or 2 indicates that 
extra-ecliptic target conditions are desired, in which the 
iterator dependent variable triggers Y1 (2) through Y6{2) are 
set equal ro 1, and for which the input LAUNCH (which see) 
should probably be set to 1, and parameters related to LAUNCH 
also set appropriately. Ordinarily M0PT2 = 3. No retro stage 
may be employed. 

=1 i, e, r p specified; f n - 0 . 

=2 i, e, a specified; f optimized. 

In the above, i = final extra-ecliptic inclination, e = final 
eccentricity, r p = final perihelion distance, a = final semi- 
major axis, and f = true anomaly at the final time. Final fl 
and a) are optimized in both cases. Related to AE, AR, and AAI. 

Numerical integration option (currently not used). 

Primer-origin-proximity step-size-conti ol indicator. Value of 
zero causes the bypass of control, leaving the step-size Au 
constant. See discussion in the section. Integration 
(Thrust) of [1]. 
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I ROT 

ISPIN 

ITF 

ITPRNT 

*JA 

**JPP 


JPRINT 

JT 


A non-zero value of IROT causes the input ecliptic projection of 
the primer vector and its time derivative to be rotated about the 
z-axis through an angle equal to the difference in longitudes of 
the spacecraft between the last trajectory of the previous case 
(or zero if no previous case) and the first trajectory of the 
current case. This feature permits one to use the initial adjoint 
variables from a 2-dimensional trajectory as the initial -guess 
inputs for a 3-dimensional trajectory using the ephemeris option. 

Spinner indicator. Not used at present. 

Provides normal termination conditions for runs which require 
more machine time than is estimated. The value specifies the 
number of machine-time seconds (CPU and I/O) required to exe- 
cute the summary trajectory after halting the iteration-sequence, 
[sec] Does not apply if subroutine REMTIM is dummied. 

Indicator for special print from MINMX3 iterator. Non-zero value 
invokes print. 

Jettison indicator j a for solar arrays (or other power generation 

system) prior to primary-target retro-maneuver, as used in 
expression (5). 

= 0 Solar arrays rot jettisoned. 

= 1 Solar arrays jettisoned prior to retro-maneuver. 

Jettison indicator j^ tor electric propulsion thrust subsystem 

prior to primary-target retro-maneuver, as used in expression (5). 
In the old spacecraft model, jettison indicator j for entire 

r 

electric propulsion system prior to primary- target retro- 
maneuver, as used in expression (9) of [1] . 

= 0 Propulsion system not jettisoned. 

= 1 Propulsion system jettisoned prior to retro maneuver. 

Unit 11 printout-length indicator. A value of zero causes the 
iterator independent and dependent variables to be output only 
for each summary-trajectory; a value of one causes the same 
output additionally at each iteration of an iteration sequence. 

Jettison indicator j^ for electric propulsion tankage prior to 
primary- target retro-maneuver, as used in expression (5). 

- 0 Tankage not jettisoned. 

= 1 Tankage jettisoned prior to retro-maneuver. 
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KPART 


LAUNCH 

LOADX 


*MATCH 


MAXHAM 

MBOOST 


Option for automatically selecting improved independent parameter 
perturbations for generating the iterator's partial derivative 
matrix. The option is invoked by setting KPART = N(N>0) , where 
N is the maximum number of allowed steps, as discussed in the 
section. Perturbation Step Size Selector of [1]. KPART must be 
set back to zero if option is not desired on subsequent cases. 

Launch mode selector, pertaining to the optimization of the de- 
parture asymptote declination, invoked by LAUNCH = 1. Related 
to X10, Y10, XI 7, and Y17. 

Intermediate-target initial-guess feature. Should be used with 
NSET (5) = 1, and then set to zero on the subsequent case. A 
non-zero value of LOADX will. invoke this feature, whereby the 
primer A and its derivative A will be loaded into the iterator 
independent- variable arrays it each intermediate-target provided 
that the trigger of the independent variable is on. The sole 
purpose of this capability is merely to generate an initial- 
guess for a multi pie- target mission, where the values loaded 
into the iterator arrays represent continuous A and A at each 
target. 

Logical indicator for maximum useable solar array (or other 
power source) power output under reference conditions (normal 
to sun at 1 AU distance) perfectly matched to maximum power 
acceptable by the power distribution system. Causes input PAO 
to be overridden. 

= T Power is matched. 

= F Power is not matched; input PAO is used. 

Maximum number of times that the program will print the warning 
message BAD HAMILTONIAN on any given computer run. 

Launch vehicle selector. 


=0 ATLAS (SLV3X)/CENTAUR 

1 TITAN III C 

2 TITAN III C (1207) 

3 TITAN III X/ CENTAUR 

4 TITAN III X (1207) 

5 TITAN III X (1207) /CENTAUR 

6 SATURN IB/LM 

7 SATURN IB/ CENTAUR 

8 SATURN IC/SIVB/CENTAUR 

9 TITAN III X (1205) /CENTAUR 

10 TITAN III B (C0RE)/CENTAUR 

11 TITAN III D (1205) /CENTAUR 

12 DELTA 

13 TITAN III D 

14 TITAN III D (1205)/ CENTAUR/TE364(2250 ) 

(continued on next page) 
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MDAY 

MODE 

MONTH 

MOPT 


MOPTX 


MOPT2 


15 TITAN III E/CENTAUR . 

16 SHUTTLE/TRANSTA6E 

17 SHUTTLE/DELTA 

18 SHUTTLE/ AGENA 

1 9 SHUTTLE/CENTAUR 

20 SHUTTLE/CENTAUR/BURNER II (2300) 

21 SHUTTLE/ I US 

MEG Use input booster coefficients B1 , B2, and B3. 

Day-of-month of reference date (e.g., 26). Related to MYEAR, 
MONTH and HOUR. 

Power variation option selector. The value of MODE is equal to 
the option-number of the power-curve, discussed in the section. 
Electric Propulsion System of reference [1] (which see). 

Possibly related to ASOL and GAMMAX. MODE = 1 has been eliminated. 

Month-of-year of reference date (e.g., 8). Related to MYEAR, 
MDAY, and HOUR. 

Ballistic option indicator. Using MOPT invokes option 2, dis- 
cussed in the section. Ballistic Trajectory Option of [1], as 
follows: 

= 0 No action (use input A . A , and v ) . 

= 1 Generate ballistic solution with flyby end conditions. 

= 2 Generate ballistic solutuon with orbiter end conditions. 
Related to REVS. 

Array of five elements, the first three of which may be currently 
used. This array specifies the target-number, or planet- 
number, of the successive intermediate- targets, and a value of 
zero indicates absence of the intermediate- target. A zero- 
entry must not precede a non-zero entry. Planet selection is the 
same as for M0PT2. MOPTX(l) pertains to iterator parameters 
X41-X50 and Y41-Y50; MOPTX (2) pertains to X51-X50 and Y51-Y60; 
and M0PTX{3) pertains to X61-X70 and Y61-Y70. Times at the 
targets are X48, X58, and X68. Not to be used unless M0PT2 / 0. 

Launch planet number and ephemeris-option indicator. 

= 0 Analytical planetary ephemeris is not used. 

f 0 Analytical planetary ephemeris is used and the specific 
launch planet is selected as follows: 

(Continued on next page) 
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=1 

Mercury 

=29 

Flora 

2 

Venus 

30 

Achilles 

3 

Earth 

31 

Amor 

4. 

Mars 

32 

Hidalgo 

5 

Jupiter 

33 

Alinda 

6 

Saturn 

34 

Grigg-Skjellerup (1977)' 

7 

Uranus 

35 

Kopff 

8 

Neptune 

36 

Grigg-Skjellerup (1982)' 

9 

Pluto 

37 

Ganymed 

10 

Ceres 

38 

Ivar 

n 

Input Target** 

39 

Bei ra 

12 

D' Arrest (1982)* 

40 

Kepi er 

13 

Encke (1980)* 

41 

Giacobini-Zinner (1985)' 

14 

Icarus (1978)* 

42 

Borrelly (1987)* 

15 

Eros 

43 

Tempel II (1988)* 

16 

Geooraphos (1983)* 

44 

Tempel II (1983)* 

17 

Encke (1977)* 

45 

Tuttle-Giacobini-Kresak 

18 

Encke (1984)* 

46 

Schaumasse 

19 

Encke (1987)* 

47 

Honda-Mrkos-Pajdusakova 

20 

Hal 1 ey 

48 

Giacobini-Zinner (1979) 

21 

Betulia 

49 

Icarus (1987)* 

22 

Toro (1983)* 

50 

Toro (1987)* 

23 

Pallas 

51 

Geographos (1987)* 

24 

Juno 

52 

Grigg-Skjellerup (1987) 

25 

Vesta 

53 

Pons-Winnecke (1989)* 

26 

Astraea 

54 

Reinmuth-1 (1988)* 

27 

Hebe 

55 

Encke (1990)* 

28 

Iris 




M0PT3 Planet number of primary target. Planet selection is the same 

as for M0PT2. If ephemeris is not used, M0PT3 is used only for 
retro-stage mass computations. 

M0PT4 Array of ten elements, specifying up to ten post-swingby targets. 

Planet selection is the same as for M0PT2, and a value of zero 
indicates the absence of a post-swingby target. A negative 
value in M0PT4(1) selects multiple ballistic swingbys, rather 
than a set of single swingbys in which case also set MAXHAM s - 0. 
Negative values (in absolute value) produce planet selection 
the same as for M0PT2 . When M0PT4(l) < 0, the remaining elements 
of M0PT4 ( i ) may be positive or negative. See the section, Swingby 
Continuation Analysis of [1] for details ana Sample Case H of [1] 
for an example-case. Should be used only for primary-target flyb.y 
missions. Related to T2, MSWING, NSWING and XSWING. 


*Year-va1ue indicates apparition for which internal orbital elements are 
most accurate. 

**Input corresponding orbit elements (see CNI, CNIX). Hone are available for 
the launch planet. 
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MPOW Flag used in conjunction with the solar array degradation option. 

Value of zero results in the optimum orientation of the arrays 
relative to the sun line. A non-zero value forces the arrays to 
an orientation yielding the maximum power achievable at that 
instant. Related to TPOWER. 

MPRINT Indicator for printing the summary-trajectory (final trajectory 
of a case) as a function of time or for invoking extra printout. 

= 0 Small-size block print at thrust switch points only 
(SWITCH POINT SUMMARY page). 

= 1 Same as = 0, except expands to become a standard print- 
block of parameters for each computed point along the tra- 
jectory, including the trajectory extension controlled by 
the input variable TGO. 

= 2 Same as = 0, except each block contains extra lines con- 
sisting of target-relative coordinates and target magnitudes. 

= 3 Combination of = 1 and = 2. 


MPUNCH Punched-card and trajectory-tape generation control. 


= 0 No special output. 

1 Punch final values of independent parameters. 

2 In addition, punch selected mission analysis parameters 
used for graphic documentation or other purposes. 


<0 and>-100 Punch trajectory output used with the ASTEA program. The 
absolute value of MPUNCH determines the frequency of tra- 
jectory points output, e.g., -3 would result in the punching 
of every third integration point. 

€-101 Trajectory tape output used with the ASTEA program. The 

absolute value less 100 determines the frequency of trajectory 
points output. Related to NTAPE. 

MREAD Card input option (iterator independent variables) 


= 0 No special cards input. 

= 1 The independent variables generated by a previous run by the 
MPUNCH = 1 or 2 option are input following the NAMELIST case, 
as discussed in the section. Program Output of [1]. 


MSWING Array of ten elements, used only when running multiple-target 
ballistic swingbys, such that MSWING(i) corresponds to M0PT4(i) 
and selects the type of swingby maneuver desired at the respective 
swingby target. Used only if M0PT4(1) < 0. The shooting method 
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(MINMX3 iterator) is used, and values of -1, -2, or -3 
correspond to a swingby passage distance initial guess of 
rp = » (i.e., continuous heliocentric velocity). Each element 

MSWING(i) may have any of the following values: 

= -1 Go* directly for unpowered swingby; if and only if it fails, 
go for powered swingby having flight time T2(i) = initial 
guess. 

= -2 Go directly for powered swingby only, having T2(i) = flight 
time of post-swingby leg. 

= -3 Go directly for unpowered swingby; then, whether it Succeeds 
or not, go for powered swingby having T2(i) = flight time. 

= -4 Go directly for unpowered swingby, but using initial velocity 
guess loaded into XSWING(j, i), j = 1,2,3, similar to 
MSWING(i) = - 1. 

= -5 Same as = -2, except use initial guess as in = -4. 

*"Go for" means "attempt to obtain (solution)". 

Related to M0PT4, T2, XSWING, and NSWING. 

MTMASS Mission-type selector pertaining to the primary target. 

- 0 Flyby mission. 

1 Orbiter (high-thrust retro-maneuver without velocity loss). 

2 Orbiter (high-thrust retro-maneuver with velocity loss). 

3 Specified arrival excess speed v mn . 

If v „ = 0, rendezvous mission 
■»n 

If v mn > 0, controlled flyby mission 
Mo retro-maneuver in either case. 

4 Orbiter (Electric propulsion system performs spiral maneuver. 
Arrival excess speed v must be specified as zero). 

Other parameters which may be related to MTMASS are DMRETR, CTRET, 
RPER, RAP, THRET, SPIRET, <JPP,JT, and JA. 

MUPDAT Flag indicating whether iterator independent variables at end of 
one case are to be updated for use as first guesses of next case. 

= 0 Do not update independent parameters. 

1 Update independent parameters for next case to be those 
obtained at end of iteration on the current case. 

MYEAR Year of reference date (e.g., 1982). Related to MONTH, MDAY, 

and HOUR. 
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*NDELTA 

NDIST 


*NEW 


NHUNG 

*NMAX 

*NMIN 

NORMAL 

NPERF 

NPRINT 


Increment An^ in number of thrusters which may be switched on 
or off at any given time. Related to NMIN and NMAX. 

Identification number of celestial body to be used as the 
reference for the communication distance and angle measurement 
printed in the Extremum Point Summary Table. Identification 
code is the same as for M0PT2. Useful for determining minimum 
distance of spacecraft to other bodies in the solar system, 
including the primary target when attempting to generate a 
solution for the first time. 

Master logical indicator for invoking the improved spacecraft 
model logic; related to all other inputs flagged by an 
asterisk(*). 

= T Use new spacecraft model 
88 F Use old spacecraft model 

Maximum number of propulsion-corner-proximity occurrences 
allowed in a given iteration-sequence. Related to GAP. 

The maximum number of operating thrusters, n . Related to 
NMIN and NDELTA. max 

The minimum number of operating thrusters, n^: greater than 

zero — pertains to thrust phases only. Related to NMAX and 
NDELTA. 

Automatic adjoint-variable scaling. 

- 0 No action. 

1 All A and A are scaled such that X „ becomes unity. 

Identification number of end condition that is to be used as 
the performance index when employing the direct parameter 
optimization feature (improve Mode). The identification code 
is the same as the i in the Yi end condition array. 

Print selection flag. Permits selection of amount of printout 
desired on each case. 

= 0 Print only the case summary. 

1 Print switching point summary of final trajectory. 

2 Print MINPUT and case setup. 

4 Print trajectory summary on each iteration. 

8 Print partial derivative matrix each iteration. 

(Continued on next page) 
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NSET 


NSWING 

NSWPAR 


NTAPE 

OMI 


Combinations of options obtained by summing options desired. 

If NPRINT >15, printout consistent with NPRINT = 0 is obtained. 
If the sign of NPRINT is reversed to negative, the iterator 
independent and dependent variables additionally are printed 
for every trajectory which HILTOP generates (including neighbor- 
ing trajectories). 

Iteration-sequence control array. 

NSET(l) Not used for input. 

NSET (2) Not used for input. 

NSET(3) Maximum number of iterations permitted in attempting 
to satisfy constraints in satisfy mode. If zero, no 
upper limit imposed. 

NSET(4) Flag indicating whether constraints are to be satisfied 
prior to entering improve mode. 

= 0 Satisfy constraints first. 

1 Proceed immediately to improve mode. 

NSET{5) Maximum number of iterations permitted after entering 
improve mode. Setting NSET(5) = 1 causes iterator to 
be bypassed and computes single trajectory to obtain 
printout. 

Swingby continuation analysis option indicator. NSWING must be 
negative and has the same definition as MSWING (which see); 
NSWING must be used when M0PT4 ( 1 ) > 0, and may be used when 
M0PT4(1) <0. If M0PT4(1) < 0 and MSWING(i) = 0, then MSWING(i) 
will be set to the value of NSWING. Related to MSWING, M0PT4 
and T2. 

Iterator independent-variable perturbation-increment control. 

= 0 No action. 

1 Allows the iterator to vary a given independent-variable 
perturbation ax whenever a neighboring trajectory is 
detected which has a different number of thrust switch 
points than the associated nominal trajectory, ax is 
varied until the same number of switch points is achieved. 

Specifies the unit-number for the ASTEA trajectory tape. Per- 
tains to when MPUNCH s - 101. 

Ascending node angle (with respect to vernal equinox) of primary' 
target orbit. Input only when M0PT3 = 11. Related to CNI, ECI, 
SAI, SOI, TPI, EMUODD, and RADODD. [deg] 

SSSttSS 


45 


OMIX 

*ONHDOT 

*PAO 

*PDIU 

*PLOSS 

*PLV 

*PMM 

*PNFLOW 

POWFIX 

*POWMAX 

*PRAMP 

*PROOT 


Array of five elements, the first three of which may be currently 
used. Ascending node angles of intermediate-target orbits. 

Input OMIX(i) only when MOPTX(i) = 11. Related to CNIX, ECIX, 

SAIX, SOIX, TPIX, EMUODX, and RADODX. [deg] 

Logical indicator which allows derivatives of hamiltonian 
switch functions to be monitored along each trajectory. The 
user should ignore this input. 

Maximum power output by solar arrays (or other power source) 
under reference conditions, p . as used in expression (2). 

aU 

Related to MATCH, [kw] 

Digital interface unit power requirement, p^. u , as used in 
expression (9). [kw] 

Thruster neutral propellant loss, m , as used in expression (14). 
[equivalent amps] 

Power processor low voltage input power requirement, p , as used 

X»V 

in expres:ion (9). [kw] 

Mission module power requirement, p mm , as used in expression (9). 
[kw] m 

Neutralizer propellant flow rate, m , as used in expression (14). 
[equivalent amps] 

Launch-vehicle-independent (i.e., no launch vehicle) trajectory 
option in which the value of POWFIX is the spacecraft's reference 
power, [kw] Not available when NEW = T. 

Maximum input power to an individual thruster, p s , as used in 

lHClA 

expression (25a). When non-zero, this input overrides CURMAX. [kw] 

Logical indicator used for special printout during implementation 
and debugging of the new spacecraft model. The user should 
ignore this input. 

General logical indicator allowing more printout from the HILTOP 
trajectory function monitoring module. 

= F No printout from the function monitoring module. 


= T Generates summary printout from the function monitoring 
module along the summary trajectory of each case. 


PSIGN 


*PTS 

RADODD 

RADODX 

RAP 

REVS 

RPER 

SAI 

SAIX 

SOI 

sorx 

SPIRET 

STATE 

STEP! 


Flag defining the sense of the launch hyperbolic excess velocity 
relative to the initial primer vector. A value of +1. results 
in the assignment of the geocentric right ascension of the excess 
velocity equal to that of the initial primer vector. A value of 
-1 . causes the geocentric right ascension of the excess velocity 
to be 180 degrees from that of the initial primer. 

Thrust subsystem power requirement, p^, as used in expression 
(9). [kw] ts 

Radius of primary target. Input only when M0PT3 = 11. Related 
to CNI, ECI, OMI, SAI, SOI, TPI, and EMUODD. [meters] 

Array of five elements pertaining to the radii of intermediate 
targets. These inputs are not used at present. 

Apoapse distance of capture orbit about primary target. 

[planet radii] 

Number of complete revolutions of the ballistic trajectory 
generated when the associated input MOPT is used. Must be a 
positive whole number. 

Periapse distance of capture orbit about primary target. 

[planet radii] 

Semi-major axis of primary-target orbit (must be positive). 

Input only v/hen M0PT3 = 11. Related to CNI, ECI, OMI, SOI, 

TPI, EMUODD, and RADODD. [AU] 


Array of five elements, the first three of which may be currently 
used. Semi-major axes of intermediate-target orbits (must be 
positive). Input SAIX(i) only when MOPTX(i) * 11. Related to 
CNIX, ECIX, OMIX, SOIX, TPIX, EMUODX, and RADODX. [AU] 

Argument of perihelion of primary-target orbit. Input only 
when M0PT3 = 11. Related to CNI, ECI, OMI, SAI, TPI, EMUODD, 
and RADODD. [deg] 


Array of five elements, the first three of which may be currently 
used. Arguments of perihelion of intermediate- target orbits. 
Input SOIX(i) only when MOPTX(i) = 11. Related to CNIX, ECIX, 
OMIX, SAIX, TPIX, EMUODX, and RADODX. [deg] 


Retro-stage specific impulse (pertaining to the retro-maneuver 
at the primary target), [sec] 

Array of six elements containing the Cartesian position and 
velocity components of the primary target. Use only when 
M0PT2 = 0 and the trigger settings of Yl(2) through Y6(2) are 
0 or 1. [AU, AU/taiO (tau = 58.132440991 days) 

Thrust-phase computation step size, au. Related to AN. 
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STEP2 

TCOAST 

TDV 


T60 


THRET 

TOFF 

TPI 

TPIX 

TPOWER 


Coast-phase computation step size, Ag . 

Array of twenty elements, consisting of the durations of the 
coast phases corresponding to the coast-phase start-times input 
in the associated array TOFF, [days] 

Time of occurrence of an impulsive deep space burn, in days 
from the start of the trajectory, which may be used only if 
the entire trajectory is ballistic (i.e., electric propulsion 
is not permitted with this option, nor is a third intermediate 
target). Iterator independent variables X64, X65, and X66 must 
be turned on, as these are used as the A V vector components of 
the deep space burn in EMQS. Also, set MAXHAM = 0. The following 
special feature is available regarding a first intermediate- 
target. If 1.D5 < TDV < 2.D5, then the burn occurs (TDV - 1.D5) 
days after passage of that target; if TDV > 2.D5, the burn occurs 
(TDV - 2.D5) days before passage of that target, [days] 

Ballistic trajectory-extension print option. When zero, no 
action. When positive, TGO = the number of days that the 
trajectory is to extend ball is tically beyond the primary-target 
when no swingby-continuation is requested, and ball is tically 
beyond the (last) post-swingby target when swingby-continuation 
is requested (in addition to the post-swingby trajectory segment 
itself). Any negative value will invoke printout of only the 
post-swingby trajectory segment or segments when swingby- 
continuation is requested. Applies also to trajectories with 
multiple swingbys. [days] 

Retro-stage thrust, f r , used only when MTMASS =2. [lbs] 

Array of twenty elements, consisting of the times, in days from 
the start of the trajectory, at which imposed coast phases are 
to begin. Times must be in ascending order. Related to TCOAST. 
[days] First negative value indicates end of input. 

Time from reference date (MYEAR, etc.) to perihelion passage, 
for the primary target. Input only when M0PT3 - 11. Related 
to CNI, ECI, OMI, SAI, SOI, EMUODD, and RADODD. [oays] 

Array of five elements, the first three of which may be currently 
used. Times from reference date (MYEAR, etc.) to perihelion 
passages, for the intermediate targets. Input TPIX(i) only when 
MUPTX(i) = 11. Related to CNIX, ECIX, OMIX, SAIX, SOIX, EMUODX, 
and RADODX. [days] 

Solar-cell degradation characteristic-time; nuclear electric pro- 
pulsion radioactive- decay characteristic-time. Related to MPOW. 
[days] The default value must be used when NEW * T. 
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*TRATIO Minimum throttling ratio, t as used in expression (25b). 

When non-zero, this input overrides CURMIN. 


TSCALE 

T2 

VOLTAGES 


XANG1 


Iterator dependent-variable tolerance- interval scaling factor; 
scales all tolerances multiplicatively by the amount TSCALE. 

Array of ten elements consisting of initial estimates of swingby- 
continuation trajectory-segment flight-times, i.e., T2(i) 
corresponds to M0PT4(i). [days] 

The beam voltages are input as iterator independent variables 
X36 through X40 (which see). 

Latitude of the launch site. Used only if LAUNCH is non-zero. 
Related to XANG2. [deg] 


XANG2 Maximum parking orbit inclination permitted by range safety 

considerations. Used only if LAUNCH is non-zero. Related to 
XANG1 . [deg] 

XSWING Array of velocity vectors consisting of initial velocity guesses 
of a given post-swingby trajectory segment. Used only when 
either NSWING or MSWING has a value of -4 or -5, See especially 
the description of MSWING = -4. Velocity consists of exactly 
the same values as found in the VI, V2, V3 locations of the 
trajectory bl _>ck print (first block). Related to MSWING, NSWING, 
M0PT4, and T2. [AU/tau] 


XO Array of seven elements, the first six of which contain the 

Cartesian position and velocity components of the launch 
planet. The seventh element is not used for input. Used only 
when M0PT2 =0. [AU, AU/tau] 


The following describes the iterator independent and dependent variable 
arrays of the boundary value problem. Input pertaining to the individual 
independent parameters is contained in the arrays XI through X70. The 
independent-parameter arrays have five elements for each variable, as 
follows (where i s 1, 2, 3, ..., 70): 


Xi ( 1 ) Input value of parameter. Must be input regardless of trigger 

setting. If trigger is on (i.e., Xi(2) = 1), input value is 
■ used as initial guess of independent parameter and is varied at 
each subsequent iteration. If trigger is off, the parameter is 
not used as an independent parameter and is not changed. 

Xi ( 2) Trigger indicating whether parameter is to be an independent 

parameter in boundary value problem. 

Xi (2) = 0 Not an independent parameter. 

(Trigger is “off"). 

1 Use as independent parameter. 

(Trigger is "on"). 
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Xi (3) Maximum change to parameter permitted in a single iteration. 

Should be a positive quantity. Used only if trigger is on. 

Units are same as that of the parameter. 

Xi(4) Perturbation increment used to compute partial derivatives by 

finite differences. Used only if trigger is on. Units are 
same as that of the parameter. 

Xi(5) Weighting factor. Should be a positive quantity. A value of 1. 

is generally recommended. The larger the weighting factor, the 
more the parameter is inhibited from varying. -'Used only if 
trigger is on. 

The independent variables are as follows: 


XI 

i 


X2 

V 2 ) 

> Initial primer vector. 

X3 

A 0 (3) J 

1 

X4 

V 1 ’ ] 

| 

X5 

V 2 ) 

> Initial primer derivative. 

X6 

h a (3) , 

> 

X7 

\o 

Initial mass-ratio adjoint-variable 

X8 

\ 

Propulsion-time adjoint-variable. Should be 
zero when NEW = T. 

X9 


Not used. 

XI 0 

6 

Geocentric declination of launch hyperbolic excess 
velocity, [deg] 


There is no conversion from input to internal units for any of the adjoint 
variables. 

2 

XI 1 Reference thrust acceleration, g. [m/sec ] 

Ignored when NEW = T. 

XI 2 Electric propulsion system jet exhaust speed, c. 

[m/sec] Ignored when NEW = T. 

XI 3 Launch hyperbolic excess speed, v m0 . [m/sec] 

XI 4 Hyperbolic excess speed at primary target, v , 

[m/sec 1 

XI 5 Initial time, t 0 , measured from the reference 

date (MYEAR, etc.), [days] 
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XI 6 
XI 7 

XI 8 
X 1 9 
X20 


Time at the primary target, t n , measured from 
the reference date (MYEAR, etc.), {days] 

Launch parking orbit inclination, i. Used only 
if LAUNCH = 1. Optimized internally by the program 
if both XI 7 and Y17 triggers are off. [deg] 

Initial spacecraft heliocentric velocity. Not 
required unless one of the three triggers is on. 
[AU/tau] (tau = 58.132440991 days) 


X21 Constant thrust cone-angle, $ . Non-zero value 

invokes the constant-^ constraint. 0 < <$ $ 180°. 
Zero-value implies that ^ is optimized aiong the 
trajectory (variable 4>). [deg] 


X22 through X29 are currently not used (although some locations following 
X21 are reserved for additional constant thrust cone-angles). 


X30 X Degradation-time adjoint-variable. Should be 

zerc when NEW = T. 

X31 through X35 are currently not used. 


X36 through X40 are the constant (along the trajectory) beam voltages, in 
volts. When NEW = T, at least X36 must be greater than zero. 


*X36 \ 

*X37 ) 

*X38 ) Beam voltages, Vj r volts] 

*X39 \ 

*X40 / 

Values of voltages must be ascending, starting with X36 as the lowest value. 
Value of zero indicates end of inputs, e.g., X39 = 0.D0 indicates there are 
three voltage levels to be simulated, X36, X37, and X38. X41 through X50 

pertain to the first intermediate target X51 through X60 pertain to the 
second intermediate target, and X61 through X70 pertain to the third inter- 
mediate target. The corresponding intermediate-target parameters are ignored 




0? so® 


5 " 


if the intermediate target is absent. Subscripts 1, 2, and 3 pertain 
to the first, second, and third intermediate targets, respectively. 
Intermediate targets are invo ked via the MOPTX array 


X41 

Vo n 

1 

X42 

V 2 > 

( Primer vector (at start of trajectory segment) 

i 

X43 

A, (3) ' 

) 

X44 

VO ] 

| 

X45 

V 2 > 

) Primer derivative (at start of trajectory segment) 

X46 

A, (3) ] 

\ 

X47 


Encounter speed at first intermediate target, 
Vi . [m/sec] 

X48 


Time at the first intermediate target, t-j , measured 
from the reference date (MYEAR, etc.), [days] 

X49 


Sample-mass factor, k sam p -j# for sample-retrieval at 
first intermediate target. 

X50 


Drop-mass factor, y for instrument-package 

dropoff at first intermediate target. 


The independent variables X51 through X60 and X61 through X70 are 
identical to X41 through X50 except that they pertain to the second end third 
intermediate targets, respectively. A third intermediate target may not be 
present when simulating ballistic missions having a deep space burn (See TDV), 
in which case X64, X65, and X66 are used as follows: 


X64 ax \ 

X65 Ay \ Deep-space velocity-increment. [Al)/tau] 

X66 az ! 
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Inputs pertaining to the individual dependent parameters are contained 
in the arrays Y1 through Y70. The dependent-parameter arrays have three 
elements for each variable, as follows (where i = 1, 2, 3, 70): 

Yi(l) Desired value of the dependent parameter. 

Yi (2 ) Trigger. If off (i.e., equal to zero}, the para- 

meter is ignored and is not considered a dependent 
parameter. Then the other two inputs pertaining 
to the parameter need not be input. If trigger is 
on, (i.e., not equal to zero), the parameter is 
considered to be a dependent parameter or con- 
straint. Certain of the parameters may have up 
to three non-zero trigger settings. These will be 
discussed individually below. 

YT (3) Tolerance of desired value (full interval width). 


It should be noted that the transversal ity conditions, which comprise some of 
the parameters, are developed under the assumption that all constraints are of 
the point constraint type. Therefore, the satisfv-mode is sufficient in 
solving any optimization problems for which a complete set of transversal ity 
conditions is available. 


The dependent-parameter arrays are as given below. T(x) represents 
"the transversal ity condition associated with x" and the function T(x) 
will have different values depending upon the constraints imposed on the 
problem. See NOMENCLATURE for definition of symbols and subscripts. 



Trigger 1 


Trigger 2 

Yl 

Ax n IAU] | 

a[AU] 

Solar distance*[AU] 

Y2 

AY n [AU] | 

e 


Y3 

AzJAU] | 

i [ deg ] 


Y4 

AX n [AU/tau) 

T(O) 

T <*„> ) 

Y5 

AyjAU/tau] 

T(«) 

T(y ) > optimal 




. flyby 

Y6 

Az n [AU/tau] ' 

T(f) 

T(i„)J 


*Appl i cable only for two-dimensional motion in the xy pi 
M0PT2 = 0* 



UJ 


ane. Also requires that 
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Under Trigger 1 above, the first set of conditions applies to ordinary 
targeting conditions for position and velocity, and also to extra-ecliptic 
conditions to be satisfied when IOUT = 1; the second set of conditions 
applies to extra-ecliptic missions when IOUT = 2. T(n), T(io), and T(f) 

are symbols for the transversal ity conditions yielding optimum final node 
angle, argument of perihelion, and true anomaly, respectively. 



Trigger 1 

Trigger 2 

Trigger 3 

Y7 

v n 

Sm 

m net [k S ] 

Y8 

T(t) 

x[days] 

Not used when NEW = T. 

Y9 

Currently not used. 


Y10 

T(fi) 

6[deg] 

Used only if LAUNCH f 0. 

Yll 

T(g) 

0 

g[m/sec ] 

P re f[kw]- Not used when NEW = T. 

Y12 

T{c) 

c[m/sec] 

Not used when NEW = T. 

Y13 

T(V„ 0 ) 

v «o [m/sec] 


Y14 

T < v .n> 

v »n tni/sec] 

extra-ecliptic inclination [deg] 

Y15 

T(t 0 ) 

t Q [days3 


Y16 

T(t„) 

t n Edays] 

t n -t Q [days]* 

*Time transversal ity with flight time fixed 
Trigger 1 . 

is assigned to Y15 under 

Y17 

T{i) 

i [deg] , 

where i = parking 
orbit inclination. 

Used only if LAUNCH f 0. 

Y18 

T(X 0 ) 

x Q [AU/tau] 


Y19 

T(y 0 ) 

y Q [AU/tau] 


Y20 

T(z 0 ) 

z 0 [AU/tau] 


Y21 

TU) 

[deg] for <(> = 
constant with time. 



Y22 through Y29 are currently not used. 
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Y30 


T(s) 


s[days] 


{Degradation time, 
when NEW = T.) 


not used 


Y31 through Y35 are currently not used. 

Y36 through Y40 pertain to the set of constant beam voltages: 


*Y36 

*Y37 

*Y38 

*Y39 

*Y40 


Trigger 1 


T(Vj) 

(Caution: see discussion 
in section 2.2 before 
attempting to optimize 
the voltages) 


Trigger 2 


Vj [volts] 


Y41 through Y50 pertain to the first intermediate target: 



Trigger 1 

Trigger 2 

Y41 

AX-j [All] 


Y42 

Ay-j [All] 


Y43 

Az-j [AU] 


Y44 

Ax-j[AU/tau] ' 

T <V \ 

Y45 

4y-| [AU/tauj | 

T(ip 

Y46 

Az^AU/tau] 1 

T(z, ) ) 

Y47 


V wl [m/sec 1 

Y48 

ntj) 

todays] 

Y49 

m samp l [k 9 ] 


Y50 

m drop l Ik9] 



optimal 

flyby 


Y51 through Y60 and Y61 through Y70 are identical to Y41 through Y50 except 
that they pertain to -the second and third intermediate targets, respectively. 
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3.3 Default Values of Input Parameters 


The following is a complete, alphabetical list of the default values of 
program input quantities having non-zero (and non-false) default values, 
except for the iterator arrays. All other inputs are zeroed (or set false). 
The default values of the iterator arrays Xt(l), Xi(2), Yi(l), and Yi(2), 
V.for i - 1,2, 3,. ..,70, are zero, and the default values of Xi (3) through Xi (5) 
and Yi (3) for the same range of i are listed in the listing of program 
inputs of Sample Case A. Exceptions to the setting of Xi (1 ) to zero are 
displayed below. 


ALPHAA 

15. 

NDELTA 

1 

ALPHAT 

15. 

NDIST 

3 

AN 

1.5 

WUNG 

25 

AR 

1 . 

WAX 

4 

BI 

.76 

NMIN 

1 

CSEP(1] 

.987 

NPRINT 

7 

CSEP(2] 

-0.018 

NSET(3) 

300 

CSEP(3) 

.08 

NSET(5) 

300 

CSEP(4) 

2.3 

NSWPAR 

1 

CSEP(5J 

.002 

NTAPE 

17 

CSEP(6) 

.9316667 

PAO 

35. 

CSEP(7) 

.0033333 

PDIU 

.02 

CTANK 

.03 

PLOSS 

.267 

CTRET 

1/9 

PLV 

.03 

CURMAX 

2. 

PMM 

.4 

CURMIN 

.5 

PNFLOW 

.04 

CVOLT 

10. 

POWFIX 

- 1 . 

DEFFIC 

.89 

PSIGN 

1 . 

DI 

13. 

PTS 

.2 

DLOSS 

198. 

RADODD 

1 . 

DVOLT 

36. 

RAP 

38. 

ETAPD 

.95 

RPER 

2. 

ETCONV 

.95 

SAI 

1 . 

GAMMAX 

1 . 

SPIRET 

300. 

GAP 

.0001 

STATE (1) 

1 . 

HOUR 

12. 

STATE(5) 

1 . 

IRK 

1 

STEP1 

.03125 

IRL 

1 

STEP2 

.125 

ITF 

3 

TDV 

- 1 . 

MAXHAM 

5 

TGO 

- 1 . 

MDAY 

1 

THRET 

400. 

MODE 

4 

TOFF 

20*-l . 

MONTH 

1 

TPOWER 

10.**3l 

M0PT3 

10 

TSCALE 

1 . 

MUPDAT 

1 

T2(i) 

50*i 

MYEAR 

1975 

X0(1) 

(Continued 

1 . 

on next page) 
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X0(5) 

1.0000015 

XI 1 

4.xl0~ 4 

XI 2 

3.xl0 4 

X36 

1600. 

X37 

2000. 

X38 

2400. 

X39 

2800. 

X40 

3200. 


SSSSS* 
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4.0 SAMPLE PROBLEM 


The sample problem consists of a Saturn orbiter mission simulation. 
The new HILTOP spacecraft model is employed. The assumptions, solution 
characteristics, and a graph depicting the optimal thruster subsystem 
operating time-history are displayed on the following pages, followed 
by a display of the complete program output. 

The complete NAMELIST input data set used to generate this sample 
problem is reproduced below. 


&M INPUT XI ( 2 ) =1 . DO , X2 ( 2 ) = 1 . DO , X3 ( 2 )=1 . DO , X4 ( 2 )=1 . T)0 , X5 ( 2 )=1 . DO , X7= 1 . DO 
X6(3)=l.D0,X13(2)=l,P0,X14(2)=l.D0,yi(2)=l.D0,Y2(3)=l.P0,Y3(2)=l.D0 
Y4(2 1=1.00 , Y5(2)=1.D0 , YG(2 )=l,BO ,713(2)=! . DO , Y14(2)=1.D0 
HOPT2=3 , MOPT3=6 , RPER=1 . IDO , RAP=11 . DO , THRET=6 . D1 , ,TPP=1 , ,TT=1 , JA=1 
MYEAP=1985 , MONTF^P , HnAY=25 , X15=82 . DO , X1G=1 P82 . DO , MBOOfST=21 , FTMASS=2 
KEW=T , NMAX= 3 , MATCII=T, ALPHAT=22 . 5D0 , PLV^O . DO , PMH=0 . DO , PDIIT=0 . DO , PTS= 0 . DO 
X1--9 .423752245854D 01, X2= 5. 55974044 575BD 00, X3= 2, 384470211509D-01 
X4=-l. 5257R1232R11D 01, X5=-7. 71564669:1 990D 01, X6= 1.0818234571050-01 
XI 3= 4.650P35939200D 03, XI 4= 5. 8140 10834 324D 03 fcEND 


rmowm page blank 
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SAMPLE PROBLEM 


1900 Day Saturn Qrbiter 

Assumptions 

• Launch on 16 December 1985 

• Shuttle/IUS launch vehicle 

• l.lxll Saturn radii capture orbit 

• Retro stage with 300 seconds specific impulse and 60 pounds 
thrust; include velocity losses in calculations 

“ Allowable beam current range, 0. 5-2.0 amps 

• Permit 5 discrete values of beam voltage - 1600, 2000, 2400, 
2800, and 3200 volts 

• Assume maximum of 8 operating thrusters, minimum of 1 thruster, 
may be switched in increments of 1 thruster 

• Solar array sized to accommodate 8 thrusters operating at 
maximum beam current and beam voltage at 1 AU 

• = 15 kg/kw; = 22.5 kg/kw; k t = .03, = 0 , 

m rs = 0> k rt = * in 

• Entire SEP system jettisoned prior to capture orbit insertion 

• C 1 = .987, C 2 = - .018, C 3 = .08, C 4 = 2.3, C g = .002, 

C 6 = .9316667, C 7 = .0033333 

Vq = 10, nQ — .89, Gj - 198, AV| - 36, np^ - .95, ~ «95 

• = .267, A n = .04, p diu - P tv - P ra - P ts » 0 

• Departure and arrival excess speeds are optimized. 


SAMPLE PROBLEM 


1900 Day Saturn Orbi ter 


Solution Characteristics 


m o 


4237 

kg 

m a 

= 

923 

kg 

m th 

= 

1384 

kg 

m p 

= 

853 

kg 

m t 

:= 

26 

kg 

m r 

= 

655 

kg 

m net 

- 

396 

kg 

v -o 


4651 

m/sec 

v »n 

= 

5814 

m/sec 


total burn time = 778 days 
average thruster operating time = 7359 hours 
specific impulse = 4978 seconds @ I B = 2 amps, V j = 
overall efficiency = .6965 * 2 amps, Vj = 3200 


3200 volts 
volts 
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• 0*0 

• 0.0 

* 

0.0 



0.0 

,xo* 

1. 0000000 00 oooooo 

, 0*0 

, O.o 

, 

0*0 



1 ."OOCCIBOCOOOOOO 

, 

0*0 


■ 0*0 

, IBAL- 

O.INTPR= 


o,icxjt= 


OflRKa 1, 

I HL= 1, 

i i © h3 

1PGT* 

0, [SPIN 

a 0 , tTFi= 

3,1 TPRNT^ 

0, JA= 

1 *JPP= 


1, JPR1NT= 

0 , JT — 

£ > 


1 ,KPART = 

0 .LAUNCH “ 

O,L0ADX= 

O.MAXHAN 

1 = 

5, M0OQST* 


21,HDAY= 

2St MODE* 

fe Q 


*♦ ,HCNTM = 

9. HCPT= 

0 .HOPTX* 

Ot 

0, 

0, 


0, 0.MOPT2* 3, 

1! C & 

HOPT 3= 

6 , HOPT 4= 0, 

0* 0, 

0# 


0, 

0, 

0, 

0* 0, 



0 ,MFOW= 

0 • MFRlNT= 

0,HPUNCH= 

0 *MR£AD= 

0 *MSWING= 


0, 0, 

Op 

' *1 Er 


0 i 

Of 0 * 

0* 0, 

0, 


0, HTMASSa 


2,MUP0AT= 

1 pMYEAR* 



I905*NDELTA= 

l,NDISTi= 

3.NHUNG= 

25,NHAX^ 


0, NM | N - 


i , NORMALS 

0, MPERFir 



0 . HP ft 1 NT* 

7,NSCT= 

0, 0. 

300 , 


0, 300# 

N5WJNG= 0,NSWPAR* 



1 *NT APE* 

1 7,FATCK=T,NEW=T 

,QNHDOT*F*PRt' v;P=P 

.PROOTaF 







*EMJ 


ORIGINAL PAGE) I& 
OP POOR QUALITY 


■ i TTiwrirTPi i) ii iii ' m ii 


CASE 1 


ITERATOR 


INDEPENDENT VARIABLES 


NC* INDEX 


VALUf 


STEP LIMIT 


PErU JR BAT l ON 


WEIGHT 


1 1 

2 2 

3 3 

A 4 

5 5 

6 6 

7 13 

a 14 


— 9*423 75224505 4 C 000 D 01 
5 *55974044575090000 00 
2*38447021150900000-01 
*1*52575123251100000 01 
—7 *71564669199< ir| 00 r >0 01 
1 .08152345716500000-01 
4.65093553919999900 03 
5.51401053432400000 03 


3.00000000000000GCO OC 

3.00000000000000000 o.i 
3 .ooooocooooodooooj OO 

3.00000000000000000 00 

3.oroooooooooocoooD on 

3.0000000000000QDOD 00 
5*00000000000000000 02 
S.OOOOOO^OOOOCOQOOD 02 


1-00 0000 0 00 OOOCO 00 t>- 0 B 

1 .oroooooo 3005 0 3 00 [>-05 
1 * 00000000000000000-00 
1 *00000000 0033 00 03 D-OB 
,1.00 000 OOOOODOOOOOCHCB 

1.00 0000000000 0000 0-05 
- 9*99999999993999 03 0-05 
9*9999*99999999993 0-04 


t .OOOOOOCOOOOCOOOOO 00 
l *r ocooor-k^ro POfmao oo 
l *00000000000000000 00 
1 *00300000000000000 00 

i -oooooooonoooooooo 00 

1*00000000000000003 00 
1 .00000000000000003 00 
1 .OOOOQOOOOOOOOOOOO 00 


DEPENDENT VARIABLES 


NO. INDEX 


VALUE 


TOLERANCE 


cn 

CTi 


1 1 o.o 

2 2 0.0 

3 3 0.0 

4 4 0.0 

5 5 0*0 

6 6 0.0 

7 13 0.0 

5 14 0*0 


9* 99999999959999900-0 £ 
9* 99999999999999900-05 
9.99999999995999900-05 
9 • 9 99 9999999 9V9y 9 00—0 5 
9, 999999999999 99900-05 
9* 999 9 9 959 9 9 999V 900— 0 ST 
9 u 99999yV99V9V 99900-0 5 
9* 99999999999999900-05 


( PflE— I TERAT ION) 


MODULAR SPACECRAFT PARAHETtRS 


5 DISCRETE VOLTAGES (VOUTSI = l OQO v 0 2000*0 2400*0 2600*0 3230*0 

CCNPUTLD MAXIMUM INPUT POWER TO AN INDIVIDUAL THRUSTER = 6*6257 KW 

CCVPUTED MINIMUM THROTTLING RATIO ^ 0-2500000 
UNIT THRUSTER REFERENCE FCWEft = ’'•66S K W 

MAXIMUM USABLE ARRAY POWER OUTPUT = 61*512 KW * 100.00 PERCENT OF AVAILABLE POWER £ MATCHED) 

HOUSEKEEPING FGWER tKw) MINIMUM = 0*0 MAXIMUM * 0.0 * INCREMENT * 0*0 

REFERENCE SPECIFIC IMPULSE “ 4S7EU4 SEC 


THIS CASH IE CONVERGED* 

I TRAJECTORIES WITHOUT PARTIAL DERIVATIVES AND 0 TRAJECTORIES WITH PARTIAL DERIVATIVES REQUIRED FOR THIS CASE* 


OISE 


I 


SWITCH PCINT SUMMARY 


PAGE 


TIME 

SEMI-MAJOR 

AXIS ECCENTRICITY 

INCLINATION 

node 

AUG PCS 

RMAG 

TRAVEL 

R1 

R2 

R3 

VI 

V2 

V3 

HASS RATIO 

THRUST ACC 

LI 

L2 

L3 

L4 

LS 

L6 

L7 

HAM 

LG 

LC 

LPH1 

CONE 

CLOCK 

HU AG 

POWER FNCT 

SWITCH FNC1 

PS I 

theta 

PHI 

LATITUDE 

LONGITUDE 

FLT PTH ANGLE 

YHAG 

P3GP TIRE 

AO* THR 

voltage 

CURRENT 

PHAN CUR 

SP IMP 

T1RUST 

BEAM POWER 

DUMP POWER 

UTIL EF 

THR KEO 

TFR EF 

PPU EF 

PTH EF 

55 EF 

ARRAY POWER 

PPU POWER 


EARTH 


START CIF TRAJECTORY. THRUST ON 


0.0 

9.450694970-02 
■9* 4 23 75 22 SO QI 
0*0 

t*254S46640-01 
6*0000000 00 00 
9*2526591 10-01 


1*516665800 00 
9*795469950-01 
5*559740450 00 
0.0 

9*213448210 01 
3*200000000 03 
9.5IC000G0D-CI 


3*521103480-01 

0*0 

2 #384470210—0 1 
0*0 

9*213447700 01 
2.000000000 00 
7.042758670-01 


1 * 926094 43 D — 02 
-1*167704500 00 
-1.525701230 01 
1*032941980 C2 
0*0 

2*041141490 00 
9.348538690-01 


■ 6*448914^20 01 
1.0 156V2 17J-01 
-7*7l56466Jt> 01 
8.5Z49j60t>0 31 
8.448914620 01 
4*970424490 03 
9*37221 J96 0-01 


-2.544443750-14 
3*944174510-04 
1*001623 460-01 
1* 15342057D 00 
-5*396683120-01 
3.945731950-01 
9.383333000-01 


9*84*954620-01 
1*000000000 00 
1 *406000001} 00 
1*02043234 O 00 
1*17211 366D 00 
5.355581480 00 
6*27689 1 078 01 


0*0 

6*983891 180-02 
2. 02 £021460-01 
9*379160960 01 
0*0 

1.053222030 PC 
7*453809090 00 


SWITCH THRUST OFF 


7.77B2 19930 02 
3*9323901 80-01 
5*949610450 00 
0*0 

— 4.3943231 CD 00 
0-L 
0*0 


6.104951960 00 
-6.184Q5402O DC 
-1.165327240 01 
0*0 

2*342100830 01 

0.0 

9*0700 CCO 00—0 1 


7.837524100-01 
8*697229120-03 
-1*056984100 00 
0*0 

2*380603890 01 

0*0 

o.o 


3*064397330-01 
2.6709SQ9GD-OI 
3*807837150-01 
7*617342910 01 
-8.087742640-02 
6*459602050-01 
0*0 


8*160037870 01 
-2.97032^420-01 
2*093624310 00 
8.139840820 01 
-8. 63497a040 31 
0*0 
0*0 


1*920504210 02 
-2. 08C16813D-03 
-5* 766502080-02 
1*534065650 00 
5* 16C23119D Cl 
0* 0 

9. 31666700 0—01 


6.176630770 00 
7 .90 569193D-0 1 
2*055816710 01 
4*0S526026D-02 
4, P* typ 1*9330-01 
0*0 

2*494536550 00 


1. B91612640 02 

0.0 

2. 027957810-01 
1*065814 100-13 
7.770219930 02 
2*494536550 00 
0.0 


SATURN 


END OF TRAJECTORY* THRUST OFF 


1*900000000 03 
4. 8594131 40 00 
1.051701140 01 
0*0 

-3*424770020 00 
0*0 
0*0 


6*104951960 00 
>8.769032440 OO 
2*422279230 01 
0.0 

1.275139740 02 
0*0 

9. 87COOOOOO— Cl 


7.8375241 00-01 
4*106390640-02 
i .626857800 00 
0*0 

1.274354590 02 
0*0 
0*0 


3.864397330—01 
1 *875501820-01 
3*508513250-03 
7.652871700 Cl 
•2.346800680-01 
4*40120 7 930- 01 
0*0 


0*160007370 01 
—2*265755080-02 
1.91 B 800200 30 
0*176695450 01 
-6.100670250 01 
0*0 
0*0 


2* 173938470 02 
-1 .273732920-03 
-8. 309283710-03 
1 . 534565 6SD 00 
3*560239400 01 
0.0 

9* 31 6667 00 D- 01 


1*002554300 01 
7*905691930-01 
2.055G1671D 01 
1.4976 87 77D— 02 
1 *889101370-01 
0.0 

9*212589390-01 


2. 14 50 46 0 ID 02 

o.o 

2*027957810-01 
5.201841620 "9 
7.770219530 02 
9*21 25 8939 D-Q 1 
0*0 


ORIGINAL PAGE & 
OF POOR QUALTTS 


CASE I 


summary 


independent parameters 


1.PIUMU-?* A237S22D 

01 1 

2,PfitM2( 

5 . 559740* D 

00 J 

3.PRIM3I 

2.3Q447023-01 J 

4.P3 0TII-1 .525761 2D 

Oil 

5.P0DT 2t — 7* 71564 C7D 

01 J 

6*Puor3t 

1*061 62350 “0 1 1 

7.LKASSC 

1 *00 00 000 D 

00 J 

0* LTAUf 

0.0 

1 

9. t 

0-0 

1 

10.DCCLN1 

0 .0 


i UACCtU 

4*0 0000000-04 ) 

' 2* V JET( 

3 . 0000 GOOD 

041 

I3.VINF0C 

4.6S09359D 

031 

14. VI NFNt 

5.01401060 

031 

15,11 ML 0 t 

B.2000000D 

01 ) 

16. T 1 HENt 

1 .90200000 

03 J 

j. 7. IPARK { 

0.0 

J 

IB.VELOI [ 

0.9 

J 

19.VEL02C 

0-0 

J 

2C.VEL03C 

0 .0 


2 1, THE X 1 ( 

0*0 

1 

22* l 

O.C 

) 

23. ( 

0*0 

1 

24. ( 

0.0 

1 

25, C 

o.a 


26, ( 

0.0 

) 

27. ( 

0-0 

1 

20* C 

o.o 

1 

29. < 

0 .0 

} 

30,UDGGR( 

0.0 


31. t 

0.0 

> 

3 2. t 

C .0 

1 

33. ( 

0*0 

} 

34, t 

O .n 

1 

35, { 

0 .0 


36 * Vb L I 1 ( 

1.60000000 

03} 

37*VCLT2{ 

2 * 0000 000 D 

031 

30. VOLT JI 

2.40000 00J 

Oil 

39. V3LT4( 

2.0QOOOOQD 

031 

4 0 * VOLT 5 1 

3 • 20000 ODD 

03 1 

4 i * PR 1 — A l 

0.0 

} 

42.PR2-A t 

0.0 

1 

43,PR3-A{ 

0.0 

1 

44, PD i-A( 

0 .0 

) 

45.PD2-AI 

0.0 


4b.P03-A{ 

0 *r 

} 

47.V1NFAC 

0 .0 

) 

4U.TIHCAI 

0*0. 

1 

49. KS AMP< 

0.0 

1 

SO.KDRUPC 

0 .0 


51.PRI-0C 

0.0 

1 

52 *PR2— D ( 

D *0 

1 

5 J ,PK3— 0 l 

0 *0 

1 

b4,P01-B[ 

0*0 

> 

5G,PD2-at 

D.O 


S6.PD3-OC 

0.0 

1 

57.V1NFC1 

0 .0 

1 

56.TIME0C 

0.0* 

1 

b9 . KS AHP( 

0.0 

1 

6H ,KO«UP( 

0,0 


6 l*Pft!“Ci 

0.0 

1 

62 . PP2-C t 

0.0 

> 

63.PR3-C t 

0.0 

1 

64*PJ l“Cf 

0.0 

J 

6S.PD2-CC 

0.0 


66* PDj-Ct 

o.c 

} 

67* V INF C ( 

0.0 

J 

06 .T1HEC { 

o.o 

} 

tiO.KSAMPC 

0*0 

1 

7C,MlRUPl 

D.O 



DEPENDENT PARAMETERS 


CTi 

VO 


l.OELTA 

X (-I.365B1D- 

OS) 

2. DELTA 

Y ( 

2. 932980- 

06) 

3, DELTA ZI 

2*443133- 

ov) 

4*D~ 

LT XD(-7.P6402D 

-09) 

5.DELT 

voc 

3.4342GD-C9 

6.DELT 

ZD 1 

1 ,364710- 

C9) 

7* 

C 

0.0 

) 

Bt 

( 

0*0 

J 

9 t 

t 

0.0 

) 

10. 

( 

0.0 

11* 

( 

0*0 

) 

12* 

C 

0*0 

1 

13, T* 

VI NFC t — 1 .240870- 

OS) 

14, T. 

VlNFNt 

“3*0i4S7D 

-OBJ 

IS, 

£ 

0.0 

£ 6* 

i 

0.0 

J 

17, 

C 

o.o 

1 

1 8 * 

t 

0 .0 

I 

19, 

{ 

0 ,r 

) 

2?. 

c 

C * 0 

21. 

l 

0.0 

) 

22* 

l 

0.0 

1 

23. 

( 

o.o 

1 

24, 

C 

0.0 

1 

25, 

( 

0.0 

26. 

C 

0.0 

) 

27* 

( 

0*0 

1 

20 o 

C 

0.0 

J 

HV, 

i 

0 .o 

) 

30, 

t 

0*0 

31, 

( 

c.o 

) 

32, 

t 

0*0 

1 

33, 

( 

0.0 

1 

34, 

i 

0.0 

) 

35, 

( 

o.a 

36, 

c 

o.o 

) 

37, 

( 

0.0 

1 

30 , 

c 

0*0 

1 

39, 

i 

n.o 

) 

40, 

t 

0 .0 

41* 

t 

0.0 

) 

42* 

{ 

0*0 

J 

43, 

i 

0*0 

1 

44, 

{ 

c *0 

) 

45, 

i 

0*0 

4 be 

i 

0.0 

1 

47. 

t 

0*0 

) 

40, 

1 

0*0 


4V, 

c 

o.*o 

J 

50, 

{ 

0*0 

61* 

i 

0*0 

J 

52, 

t 

0.0 

> 

S3* 

( 

0*0 

i 

b4. 

1 

O.C 

1 

55* 

i 

0*0 

56* 

t 

0.0 

1 

57* 

t 

0.0 

) 

bU, 

f 

0*0 

) 

59, 

c 

0*0 

} 

60* 

< 

0*0 

61, 

< 

c.o 

) 

62* 

l 

0*0 

J 

63 * 

c 

0*0 

) 

64, 

[ 

0-0 

1 

65* 

t 

0 .0 

6 6, 

( 

0.0 

) 

67, 

i 

0*0 

1 

66, 

c 

0*0 

) 

69* 

£ 

0*0 

> 

70, 

{ 

0*0 

THRUST SWITCHING TIMES 

4 DAYS 1 

0 

•0 OH 

777.022 

OFF 

1900*000 

0= 

F 








POWER 

EFFICIENCY 

ELECTRIC PKGPULSIOA 
PROP TIME 

PARAMETERS 

J 

PROP TIME RATIO 

61.5120826217 

0*656524X624 

777.6219925092 

5*1016639X14 

0*4093799961 


AVE ACCEL 
0*000 463b 053 


INITIAL 

4236*7038003202 


PROPULSION 

2306.7030903X35 


MASS COMPONENT 
PROPELLANT 
053 .*026032934 


BREAKDOWN 

TANKAGE 

25.6020004900 


structure 

0*0 


PAYLOAD 

396*4636542323 


126 THRUST COMPUTE ST EPS; 


22 COAST COMPUTE STEPS 


CAS£ 1 


EXTGfcMUH POINTS OF SELECTED FUNCTIONS 




ECLIPTIC 

SOLAN 


CCMMUUICATION 


SWITCH 


THHUST ANGLES 


INPUT 


AUUAY 

ot h rn 

I 

TIME 

LONUI rUDE 

DISTANCE 


AH (iLL 

PI 

STANCE 


FUNCT 

1 CN 

Pbl 

IHC1A 

PHI 

POWLtt 


ahull: 

TAULL 

0 

0.0 

0 .0 

O . <i 0 4 


o?,v 


0-0 

ON 

9 . JCJll 

01 

0*1 

92-1 

92.1 

62 .U 

UN 

0.0 


1 

1 * 764 

2.1 Hi N 

0 - V 6 4 


09*7 


0.00 


9-40D 

01 

0,1 

91 .5 

91 .S 

MAX 62.0 


0,0 


tz 

10.672 

21.9 

1.000 


107*0 


0*05 


9 • 4 ID 

01 

0-1 

0b. 0 

85.0 

61 .5 


0.0 

l MAX 

5 

28*131 

32.7 

1.023 


11 7.2 


0-08 


9.200 

01 

MAX 0.1 

81*3 

01*3 

59.7 


0.0 


6 

37*494 

42.9 

1*056 


1 20.1 


0.11 


9,0 TO 

01 

0.1 

77,0 

77.0 

57,3 


0.0 

- V 

6 

4 7* 0 82 

52*7 

1-053 


139.9 


0.14 


0 • 7{JD 

61 

0,1 

74*3 

74.3 

S,.3 


0.0 

I MAX 

6 

59.800 

64. S 

1*166 


155.8 


0.19 


0*350 

01 

0*1 

70.0 

70.0 

50.0 


0.0 

-V 

4 

67. SOS 

71 .3 

1*214 


165.0 


0.23 


0.0 2D 

01 

0*1 

67-5 

67.5 

47, l 


0.0 

l MAX 

0 

79-424 

60.6 

1-293 

MAX 

179.9 


0.30 


7*560 

01 

0-1 

63*9 

63.9 

42.9 


0.0 


4 

80*251 

61 m2 

1-299 


179.1 


0.31 


7 .5 30 

01 

0, 1 

63.7 

63.7 

42.6 


0.0 

-V 

4 

80- 127 

66.7 

1-355 


170-5 


0.36 


7.1 OD 

01 

0.1 

61-5 

Cl o 5 

39.9 


0.0 

I MAX 

6 

IftJ. 399 

96.3 

1-471 


155.5 


0.50 


fa *6 ID 

01 

O.l 

57- 6 

57*6 

35.2 


0.0 

-V 

5 

1 11-936 

101. 1 

1.536 


iae. 1 


0.69 


6*230 

-01 

o.l 

55*6 

55.6 

32,7 


0.0 

X MAX 

14 

126.194 

109.2 

1.669 


135*3 


0.00 


5.6GO 

01 

0.0 

52,1 

52.1 

20*7 


r .0 

-N 

4 

120*226 

109.2 

1 .669 


135.3 


0.80 


5* 7 QD 

0! 

0.0 

52.1 

52.1 

28*6 


0. 0 

1 WAX 

9 

147.600 

117.4 

1*029 


121*0 


1*08 


5*1 ID 

01 

“0 .0 

48.4 

40.4 

24.6 


0. 0 

-fa 

4 

147.777 

117. S 

1.S3C 


121*7 


1.08 


b. t 30 

r\ 

-0-5 

40.4 

40.4 

24.6 


n. ■) 

I MAX 

4 

171.965 

126.0 

2.032 


106-4 


1-SO 


4.4SD 

01 

-0.1 

44*4 

44 .4 

20.5 


0.0 

-N 

4 

172.405 

126.2 

2.036 


106*2 


1.50 


4.52D 

01 

-0*1 

44.3 

44*3 

20.5 


0.0 

I MAX 

9 

204.451 

125.3 

2.303 


U7*2 


2*12 


J.UTD 

01 

-0.2 

39*0 

39*0 

16.5 


0.0 

-N 

5 

2C5. 401 

135.6 

2.312 


06*6 


2.14 


3.0 20 

01 

-0*2 

39.6 

39.6 

l 6 , 4 


0*0 

I MAX 

0 

251 .671 

145.0 

2-690 


59-0 


3*0 0 


3.020 

Oi 

-0 .4 

34.4 

34*4 

12.6 


0.0 

-N 

4 

254*203 

146.3 

2.711 


50.1 


3.11 


3 * 0 50 

91 

-0.4 

34.1 

24.1 

12,3 


0.0 

I max 

9 

325-708 

157.5 

3.272 


13.4 


4,23 


2*1 2D 

01 

-0.7 

20.0 

20. 0 

a. 7 


n .0 

-fa* V 

e 

330.91 1 

150.2 

3.311 


10-0 


4,28 


2*270 

01 

-0.7 

27.7 

27.7 

8*5 


0*0 

-V 

4 

338. 746 

159,2 

3.370 


4-7 


4*35 


2. OSD 

01 

-0*7 

27.1 

27.1 

0.2 


0* 0 

I ha x 

5 

345.607 

160.1 

3-422 

MIN 

0-0 


4*41 


2.0 2D 

01 

-o.a 

26.7 

26.7 

6.0 


0*0 


5 

383.759 

164.4 

3.704 


27*5 

MAX 

4-50 


1 .G7D 

01 

-1*3 

24.4 

24.4 

6.8 


c* o 


6 

445*356 

170-2 

4.142 


77-5 


4.24 


1 * 1 8D 

01 

-1 .3 

21 ,6 

21 *6 

5.5 


0.0 

— N+ v 

7 

4 73*353 

172*5 

4*334 


103.6 


3.99 


1 *35D 

01 

-1 *5 

20.6 

20*7 

5,0 


0.0 


6 

476.B26 

172.7 

4.350 


10 7.1 


3,96 


1.26D 

01 

-1 . b 

20*5 

20*6 

5.0 


0.0 

I WAX 

4 

517.853 

175.7 

4*621 


15C.5 


3*72 


1 .0 ID 

Q 1 

-1*3 

19.5 

19.6 

4.4 


0.0 

+ v 

S 

523-088 

176.1 

4.665 


156-4 

MIN 

3*72 


1 *96D 

01 

-1.0 

19*4 

19.5 

4.4 


0.0 


5 

S43.992 

177.5 

4.ar i 

MAX 

100.0 


3*79 


9* 4 ID 

00 

-2.0 

19*0 

19.1 

4.1 


0.0 


0 

56 1.027 

178.5 

4.909 


160.9 


3.94 


0.5 ID 

00 

-2*1 

10 .8 

19.0 

3.9 


0.0 

*v 

9 

593*995 

100.5 

5.1 16 


126-1 


4,45 


7*3SD 

00 

-2.4 

10.7 

10.0 

3.6 


C .0 

+ v 

6 

595. 123 

180.5 

5*122 


1"4.9 


4.47 


7*760 

00 

-2*4 

10.7 MIN 

18. H 

3,6 


0.0 


6 

599,255 

180.8 

5.140 


120.0 


4.55 


7.S60 

00 

-2*4 

HIN 10.6 

16.0 

3.6 


0.0 


5 

737*074 

157.5 

5.95e 

HIN 

0-1 


6.94 


1 . 5 ID 

00 

-3.9 

21,4 

21 .7 

2.7 


n.n 


5 

760.592 

100.5 

6*083 


19.2 

MAX 

7.00 


6*3 90— D 1 

-4.2 

22.5 

22.8 

2.6 


c.o 


5 

777*022 

189.2 

6.177 


33.9 


6.97 

OFF 

1 *070-13 

-4.4 

23.4 

23.8 

2.5 


0.0 


6 

911.527 

193.9 

6-857 


163*3 

MIN 

5*80 


-4.27D 

00 

***** 

** **** 

***** 

0.0 


90*0 


5 

927.002 

194.4 

6.931 

MAX 

179.9 


5*91 


-4,600 

00 

***** 

****** 

***** 

0.0 


90.0 


7 

1 115.061 

199.7 

7,753 

MIN 

0*1 


0*74 


-8. I 10 

00 

***** 

****** 

***** 

0.0 


90*0 


6 

l 130*656 

200. 1 

7.016 


13*8 

MAX 

B.77 


-B.14D 

oo 

***** 

****** 

***** 

0-0 


90 . 0 


5 

1290.957 

203. a 

0,412 


168 cO 

MIN 

7,41 


-7 . 5 40 

00 

***** 

****** 

***** 

0.0 


90.0 


6 

1302.302 

204.0 

0*451 

MAX 

179*6 


7*43 


-7.380 

00 

***** 

****** 

***** 

0.0 


90*0 


5 

148©. 139 

207*7 

9.040 

HIN 

0.2 


10*02 


-4. 770 

00 

***** 

****** 

***** 

0.0 


90.0 


6 

1499*225 

2C7.9 

9.072 


10.0 

MAX 

10-04 


-4.54D 

00 

***** 

****** 

***** 

0-0 


90.0 


5 

1666,650 

210*0 

9.510 


171.4 

MIN 

0.51 


-7.36D-01 

***** 

****** 

***** 

0.0 


90-0 


7 

1674.021 

210.9 

9*538 

MAX 

179.8 


8*52 


-5.40D-01 

***** 

****** 

***** 

0.0 


90.0 


7 

1859.470 

£13.9 

9*946 

MI W 

0*2 


10*93 


4.120 

00 

***** 

****** 

***** 

0,0 


90,0 


6 

1 867.135 

214-0 

9.962 


7*0 

MAX 

10.94 


4*320 

00 

***** 

****** 

***** 

o.o 


90-0 


4 

1900.000 

214.5 

10.026 


37.0 


10* BO 

OFF 

S .200 

00 

***** 

****** 

***** 

9.0 

ON 

90.0 
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RO* OF 

□ E AH 


SEAM 

PM AMT UH 

UT1L1 l 

SPECIFIC 

THRUSTE? 

EXCESS 

THRUST 


T I Me 

THRUSTERS 

VOLTAGE 

CURRENT 

CURRENT 

tFFIC 

IMPULSE 

POWER 

PU* Eft 

CL 051 


0-0 

a 

32C0 


2 *0000 

2*0411 

0-9253 

49 70 

0 . 3i? 

1. 05 

0.395 


1.764 

a 

32C0 


2 .07 00 

2.0416 

0.9253 

4970 

5*35 

l. 06 

0*395 


1 3.672 

6 

3200 

* 

2 ,00 00 

2.0000 

0.9253 

4970 

b. 36 

0. 0 

0.39E 


26* 131 

6 

3200 


L .94 1 4 

1-9414 

0*9 198 

4955 

5* 18 

0. 0 

0* 383 


37*494 

6 

- 3200 


1 .0622 

1*6622 

0.9121 

4920 

T 4*94 

0. 0 

0,368 


HCPCAT 

S 

2000 


2 .0000 

2.1106 

0.9253 

4657 

4 • b ? 

2. 40 

0*305 


4 7*062 

a 

2000 

* 

2.0000 

2-0000 

0*9253 

4657 

4.6? 

0*0 

0*365 


59*600 

a 

- 2 onn 


i . a j 69 

i.ojb9 

0,9097 

4593 

4.26 

0.0 

0*340 


REPEAT 

a 

2400 


2.0000 

2*1216 

0*9253 

431 1 

4.02 

2. 34 

G> 342 


67.505 

a 

2400 

♦ 

£ .0000 

2-0000 

0.9253 

431 1 

4,32 

0, 0 

0*342 


79-424 

6 

2400 


1 .0192 

1.0192 

0.9077 

4244 

3.61 

0, 0 

0*31 £ 


00*261 

a 

- 2400 


1 .0069 

1-UC69 

0.9364 

4239 

3*58 

0* o 

0*310 


t -REPEAT 

0 

2000 


2 *0000 

2. 1347 

0*9253 

3936 

3*35 

2* 15 

0.212 


as. 127 

6 

2 CO 0 

<* 

2 .0000 

2.0000 

0,9253 

3936 

3 * 3b 

0* 0 

0.312 


103. J99 

a 

- 2000 


1 * 75 93 

1 * 7693 

0.9013 

3851 

2* 89 

0, 0 

C. 276 


REPEAT 

0 

1600 


2.0000 

2.1485 

0.9253 

3520 

2.63 

1 * 90 

0.279 


1 1 l *936 

0 

1600 

* 

2.0000 

2 .0000 

0-9253 

3520 

2.68 

0.0 

0.275 


126* 194 

- a 

itoo 


1 *7491 

1.7491 

0.9002 

3441 

2a 30 

0* 0 

C.245 


REPEAT 

7 

1600 


2.0000 

2*0005 

0,9253 

3520 

2,68 

0. 01 

0. 244 


1 26*226 

7 

1600 


2.0000 

2-OCOO 

0*9253 

3520 

2*63 

0, 0 

0*244 


147.606 

- 7 

1600 


l *7150 

1. 71b0 

C * H 963 

3429 

2.25 

0.0 

€.210 


REPEAT 

6 

1600 


2.0000 

2*0026 

0*9253 

3520 

2* 68 

0*03 

0.209 


147*777 

6 

1 too 

* 

2.0000 

2-0000 

0,9253 

3520 

2.08 

0*0 

0.205 


1 71*966 

- 6 

1690 


l ,,6702 

1.6702 

0.8912 

3412 

2.18 

0*0 

0.176 


—REPEAT 

5 

1600 


2*0300 

2.0063 

0.9253 

3520 

2.68 

0. 05 

0*174 


172.405 

5 

1600 

♦ 

2.0000 

2.0000 

0.9253 

3520 

2. 60 

0.0 

0*174 


204*451 

“ 5 

1600 


1 *6306 

1.6C66 

0.6837 

3307 

2.09 

0, 0 

0* 141 

-■J 

REPEAT 

4 

1 too 


2.0300 

2.0133 

0.9253 

3520 

2*66 

0* 09 

0. 140 

— 1 

205.401 

4 

1600 

* 

2.0000 

2.0CCD 

C.9253 

3520 

2* 68 

0. 0 

n. 140 


251*671 

- 4 

1600 


1 .5193 

1 P 61 93 

0.8723 

3349 

l* 9b 

0. 0 

0* 107 


REPEAT 

3 

1600 


2-0000 

2,0289 

0.9253 

3520 

2,68 

0, 14 

0 * tos 


254*263 

3 

1600 

* 

2. 00 CO 

2.QC0C 

0.6253 

3520 

2-63 

0. 0 

0*105 


325*768 

- 3 

♦ 1600 


1 *40 76 

1 , 4076 

0*8569 

3297 

1.78 

0* 0 

0*074 


REPEAT 

2 * 

2CC0 


1 ,7327 

1 *7327 

0*6983 

3041 

2*84 

0. o 

0*063 


33C * 9 1 1 

2 

- 2000 


1 ,6940 

1 • 694 0 

0.8939 

3025 

2.77 

D* 0 

0-066 


REPEAT 

2 

1600 


2*0000 

2.0687 

0.9253 

3520 

2.06 

0* 22 

0*070 


336*746 

2 

1600 

* 

2*0000 

2*0000 

0*9253 

3520 

2*68 

0. 0 

r,n?o 


3*5.607 

2 

1600 


1 .9420 

1.9420 

0*9200 

3534 

2*59 

0* 0 

D. Q6S 


303*759 

2 

t 600 


1 -6t?90 

1 * oo 5 0 

0*8910 

3411 

2* 1 3 

0* 0 

C* 065 


445*356 

- 2 

f* 16CQ 


1 *3439 

1.3439 

C.8474 

3264 

1*69 

o.o 

0-047 


REPEAT 

1 

2400 


1 *6701 

1.8701 ' 

0*9129 

4264 

3, 72 

0, 0 

0, 040 


473*353 

1 

- 2400 


1 .71 1 V 

1*71 11 

0.8959 

4 197 

3,3b 

0. a 

0- 037 


REPEAT 

l 

2000 


2.030' 


0*9253 

3936 

3,35 

0* 04 

0.039 

o o 

476*826 

l 

2 COO 

4 

2 ,0000 

2 *0u 0 

0*9253 

3936 

3.3b 

0, 0 

0*039 

g> 

517*053 

% 

+ 2000 


I o774Q 

1.7740 

0.9029 

3657 

2. 92 

0. Q 

0.03S 


REPEAT 

1 

2 400 


1 .50 18 

1. 5018 

0.8700 

4092 

2.09 

0. 0 

0*032 

*v yi 

523*006 

I 

2400 


1 o 46 0 1 

1.400 l 

0.8670 

4000 

2.84 

0. 0 

0-03 2 

g s 

543.992 

1 

2400 


1 .39 06 

1.3586 

0.8 656 

4032 

2.6b 

3.0 

0, 03C 

S fe 

561.027 

1 

+ 2400 


l *33 70 

1 .3578 

0.8465 

3994 

2.b2 

0 , 0 

0. 029 

£-4 

REPEAT 

1 

' 2000 


1*1599 

1 . 1599 

0.8171 

4170 

2.48 

0.0 

0*027 

JD he} 

553.995 

l 

* 2800 


1 *0607 

1,0687 

0*8000 

4097 

2 * 24 

0. 0 

0*025 

"Z_4 ' u 

REPEAT 

1 

2 200 


0.9433 

0*9433 

0.7736 

4246 

2.20 

0* 0 

0*024 

> K 

553-123 

l 

3200 


0 .9407 

0, 9407 

0-7730 

4243 

2*19 

0. 0 

0*024 

Lj-t M J 

tT* Txt 

599*255 

l 

3 250 


0.9315 

0*9315 

0*7709 

4232 

2. 16 

0 • 0 

0*023 


7J7.074 

1 

3200 


0 .6940 

0*uv48 

0.7054 

38B9 

1-49 

0-0 

0*018 

3 E? 

700*592 

1 

3200 


0 .6662 

0.6662 

0,6956 

3837 

1 , 41 

0, 0 

o,ni7 

77 7 * 822 

1 

3 200 


0 *6460 

0 * 6460 

0-6003 

3799 

1,35 

0* 0 

0*016 


91 1.S27 

0 

0 


o-o 

0-5223 

0.6371 

3524 

0*0 

0- 0 

0-0 


927-002 

0 

0 


0-0 

0.5111 

0*6317 

34 95 

o.o 

o. a 

c.o 


1115*061 

0 

0 


0-0 

0-5503 

0*6406 

30*9 

O.o 

0. 0 

0*0 


1 130-656 

0 

0 

0-0 

0 *621 6 

0.6360 

3050 

0.0 

0.0 

0.0 

1 290*957 

0 

0 

0 -0 

0.5291 

0.6403 

2799 

0*0 

0.0 

0.0 

I 302-302 

0 

0 

0.0 

0.5240 

0.6379 

2709 

0.0 

0. 0 

o.o 

14B8, 139 

0 

a 

0 .0 

0. 55b 1 

0- 662 4 

2550 

0.0 

0. 0 

0.0 

1499.225 

0 

0 

n,0 

0*5520 

0*6506 

2543 

0*0 

0.0 

0.0 

1666*660 

0 

0 

0 *0 

0.499b 

0.6261 

2450 

0-0 

0.0 

0.0 

1674.82) 

0 

0 

a .o 

0.4973 

0.6250 

244S 

0.0 

a* a 

0.0 

1 QS9.478 

0 

0 

0.0 

0.4556 

0.6031 

2362 

0*0 

0* 0 

0-0 

1667 .135 

0 

0 

0-0 

0.4641 

0.6023 

2359 

0.0 

0* 0 

0.0 

1900.000 

0 

0 

0.0 

0.4481 

0*5990 

2346 

0.0 

0.0 

0.0 


HI SSI (UN SCHEDULE 


254fisi^ii:flaa-Qa«^uuiAu w u6x£ depart earth 



X 

Y 


z 

XDOT 


vour 

ZDOT 

radius 

LAI . 

LONG * 

planet 

S/C 

9. 45 C 69500—02 
9.45C6950D-02 

9.7954 6S9D- *>! 
9.79 546990" Cl 

0.0 

0-0 


-1.01182540 

-1.1677946D 

00 

00 

9.23727930-02 
L .01569220-0 I 

0.0 

3* 9441 7450—114 

9. 84095460-01 
9.0409546D-01 

0.0 
a. o 

84 .489 
£4.409 


rv> 


24513 Jlfc*QnaQ-2Q-JULJt&H_ttaX6 PASS SATURN AT 5*8 )* KH/5EC 



X 

Y 

Z 

XDOT 

YOOT 

ZDOT 

RADIUS 

LAT, 

LONG. 

planet 

S/C 

4.U5941310 00 - 
4. £5941310 00 - 

0.76903250 

6.76903240 

00 -4*1063 9 0 90-02 
00 —4*10 639 060-02 

2.65143350-01 
1 *87550 1 60—01 

1. 56055003-01 
-2.26576510-02 

-1. 32764S2D-02 
-1.27373290-03 

1.0D2S543D 01 
1 .00255430 01 

-0.235 
— 0* 235 

-61 *007 
-61 .007 

TKO-SODY 

TRANSFER ANGLE 

BETWEEN 

EARTH AND 

SATURN 

t? 214*5049 

DEGREES. 






CASS 1 


CCCNVERGEDJ 


PERFORMANCE SUMMARY 


EARTH TO SATURN RBNDEZ WITH HIGH THRUST CAPTURE MANEUVER* 


LAUNCH VEHICLE IS SHUTTLE/IUS 


( CQEFF1 CtENTS « 209690*4200 3661.6300 3757*87971 


LD « DEC 16* 19C5, 12,0000 HOURS SKI 

JULIAN DATE 40416*0000 


AO m FEU 20* 1991* 12*0000 HOURS 3 MT 

JULIAN DATE 4U31S.0QO0 


FLIGHT TJHE * 1900-0000 DAYS, 


ALPHA A CKG/KVJ ALPHA T (KG/KVJ 

IS.0COO 22-5000 


ELECTRIC PROPULSION SYSTEM PARAMETERS 

TANKAGE FACTUR STRUCTURE FACTOR 

0*0300 0-0 


EFFICIENCY COEFFICIENTS 
0 O IKM/SbO E 

0*76000 13,00000 0,0 


ELECTRIC PROPULSION SYSTEM HASS SUMMARY { KG) 


INITIAL POWER SOURCE 

4236*7039 922*6012 


ENGINE PPOPELLANT 

1304*0219 853.4027 


TANKAGE 

25*6021 


STRUCTURE 

0-0 


NET HASS 
396*4837 


ELECTRIC PROPULSION SYSTEM PERr QUM ANC- SUMMARY 


CO 


P(1 AU) tKW) 
61*5121 


PlHSKPl IKVJ PtTARG J (KV) THRU AUl <N) ACC1 l AU1 (M/SLC + *2) 
0*0 0*9213 1 *7 SSI 49 4 *1 427ZJD~Q4 


ISP (SEC! 
4976*4 24 


EFFIC 

0*69652 


CHAR DCG { DAYS 1 
1*000 00 0 OD 30 


EXTREME TRAJECTORY AND PBRFORKANC£ CONDITIONS 


WAX OIST iAO} KIN OIST tAUI 
10 *025543 0 0*9839275 


MAX POWER tXk) 
62*702322 


MAX THRUST CN) 
t .75514909 


tiUH-4 TIME 10 AYS) 
777*8 2 199 


OEGRD TIME (DAYS) 
£48 * 50263 


TRAV ANG (DEG) 
21 4,50468 


POWER SOURCE JETTISONED PRtUR TO RtTHU MANEUVER 
ENGINE MASS JETTISONED PRIOR TU RETRO MANEUVER 
TANKAGE HASS JETTISONED PRIOR TO RETRO MANEUVER 


g ! 

nJ S 

8 | 

SO gi 


DEPARTURE AND ARRLVAL CONDITIONS 

DEP OECL ICEGJ PARK INC (OHG| OEP VINF IM/SEC) C3 tKH**2/3&C*<*2 J APR VIAF CM/SECJ 

1*4758 26*5000 4650*93694 21*631206 5814-01083 


C4 C KH**2/SEC’t 42 I 
33.802722 


high THRUST CAPTURE MANEUVER STAGE AND DTUIT SUMMARY 


STRUCTURE I KG) 
65*4512 


PERI APSE IRAQI X) 
1*1000 


PROPELLANT (KG) 
569*0611 

APCAPSL (RADII 1 

11 .0000 


THRUST (LBS) 
60*0000 

URGIT vet (M/SE C) 
33003*4812 


IS* 15ECI 
30 0.0000 

ObL VEL (M/SEC) 
£418*5247 


BURNING T (ME (SECl 
662*1310 

VEL LOSS tH/SE CJ 
322.7027 


i 



5.0 PROGRAM INSTALLATION AT LeRC 

The new version of HILTOP is stored on the UNIVAC system at Lewis 
in terms of individual source and object subroutine modules, under the 
name of MASTERS in the general disk file MFILE. A complete alphabetical 
cross-reference list between HILTOP subroutine names and source module 
names (disk file MFILE) on the UNIVAC system at Lewis is displayed 
on the next page. 

Implementation of the more sophisticated spacecraft model in the 
HILTOP computer program has enlarged the program approximately 40K bytes. 
Although the older version of HILTOP could execute within the 65K word 
capacity of the UNIVAC machine (with a rather simple overlay scheme), the 
newer version will not fit on the machine at all unless certain size- 
saving steps are first taken. In order to construct a viable overlay 
scheme for the HILTOP program, it was deemed necessary to adopt the 
ground rule that no. overlay action be allowed during a HILTOP iteration 
sequence (solution of the two-point boundary value problem). It is felt 
that if this ground rule were not adhered to, i.e., if overlay action 
occurred during an iteration sequence, the execution-time properties of 
the HILTOP program would be significantly degraded. (Indeed, it is not 
certain, at this point, if significant core savings could be obtained by 
ignoring this ground rule.) 

The following necessary size-saving steps were taken at LeRC: 


(1) A sizeable chunk of code associated with ballistic swingby continu- 
ation and ballistic trajectory extension was dummied-out. This 
consisted of writing a dummy subroutine MORE which, if referenced 
during execution, writes a simple message to the user explaining 
the non-availability of the attempted program option (s). The 
dummy M0RE is then used in the overlay structure, and this 
allows the deletion of subroutines SWING, SWTRAJ, TAPSET, and 
C0NVER from the overlay tree. It is emphasized that all dummy- 
rep! acemenrt operations have not destroyed the associated program 
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Cross reference list between HILTOP subroutine names and source module 
names (MFILE) on the UNIVAC system at the Lewis Research Center. 


HILTOP 

MFILE 

Overl ay 

Subroutine 

Module 

Segment 

Name 

Name 

Name 

AEINWT 

A 

R00T 

ALBED0 

ALBED0 

TRAJEC 

ANSTEP 

C 

R00T 

BEGIN 

D 

SBEGIN 

BHAM 

E 

TRAJEC 

CARKEP 

F 

TRAJEC 

CDERIV 

G 

TRAJEC 

CHECK 

H 

TRAJEC 

C0NVER 

I 

absent 

C0NVRT 

J 

TRAJEC 

C0RNER 

DUMC0R/K 

R00T/absent 

DATE1 

L 

R00T 

DECLIN 

M 

TRAJEC 

DERIV 

N 

TRAJEC 

EFM 

0 

R00T 

EFMPRT 

P 

FINAL 

ETA 

Q 

R00T 

ETAU 

R 

R00T 

EXTAB 

S 

R00T 

FIND 

T 

TRAJEC 

FINISH 

U 

R00T 

FUNCT 

V 

TRAJEC 

GETAMP 

W 

TRAJEC 

GETHAM 

X 

TRAJEC 

GETHOO 

Y 

TRAJEC 

GETI 

Z 

TRAJEC 

GETQ 

AA 

TRAJEC 

GETRV 

AB 

TRAJEC 

GUESS 

AC 

BALL IS 

GUNTHR 

AD 

TRAJEC 

IMPRNT 

AE 

BALL IS 

IMPULS 

AF 

BALLIS 

INC0ND 

AG 

TRAJEC 

INPUT 

AH 

INITAL 

INTERP 

AI 

TRAJEC 

L0AD 

AJ 

TRAJEC 

MAIN 

AK 

R00T 

MINMX3 

AL 

R00T 

M0RE 

DUM0RE/AM 

FINAL/absent 

NEW I NT 

NEWINT 

INITAL 


HILTOP 

MFILE 

Overl ay 

Subroutine 

Module 

Segment 

Name 

Name 

Name 

0MASS 

AN 

R00T 

PARINC 

A0 

LIMB0 

PDATE 

AP 

FINAL 

PMPINT 

AQ 

R00T 

PRINT 

AR 

R00T 

PRINTR 

AS 

R00T 

PRI0R 

AT 

TRAJEC 

PUNCH 

AU 

LIMB0 

QPRINT 

AV 

SEGQPR 

QSTART 

AW 

INITAL 

RADAR 

AX 

R00T 

RECHEC 

AY 

TRAJEC 

REMTIM 

REMTIM 

R00T 

RETINJ 

AZ 

TRAJEC 

RKSTEP 

BA 

TRAJEC 

SC0MP 

BB 

R00T 

SETUP 

BC 

INITAL 

SMQINT 

BD 

R00T 

S0LAR 

BE 

R00T 

SPRINT 

BF 

TRAJEC 

STEP 

BG 

TRAJEC 

ST0RE 

BK 

TRAJEC 

SUMMRY 

DUMSUM/BI 

R00T/absent 

SWING 

BJ 

absent 

SWST0 

BK 

TRAJEC 

SWT RAJ 

BL 

absent 

TAP 

BM 

R00T 

TAPSET 

BN 

absent 

THANGD 

B0 

TRAJEC 

TIKT0K 

BP 

R00T 

TL'iLFT 

TIMLFT 

R00T 

TRAJ 

BQ 

R00T 

TRAJI 

BR 

TRAJEC 

TRAVEL 

BS 

TRAJEC 

TWINKL 

BT 

INITAL 

VD0T 

VD0T 

R00T 

VMAG 

VMAG 

R00T 

VPRINT 

BV 

LIMB0 

VSCAL 

BW 

R00T 

WRAPUP 

BX 

TRAJEC 
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capabilities. Both the full source and object codes remain 
available for each subroutine on the UNIVAC disk pack, but the 
dummied object code is specified in the overlay tree which is input 
to the linkage editor program. The absence of the ballistic 
trajectory extension capabilities from the UNIVAC 1100 version 
of HILTOP in no way degrades the program's capability to perform 
electric propulsion performance investigations. 

(2) Subroutine C0RNER was dummied -out. This only affects simulations 
involving the "old" spacecraft model, and only when a trajectory 

is "hung" on a propulsion-time corner ( as explained in the primary 
HILTOP program document). The program will no longer automatically 
attempt to avoid this numerical difficulty by varying the propulsion 
time adjoint variable X^ , but the user can still do this manually. 
Or, when simulating the old spacecraft model, the user could use 
the older version of HILTOP which is currently being maintained 
on the UNIVAC 1100 at Lewis. Since there is no X in the new 

T 

model, this change does not degrade the program's capability when 
simulating the new spacecraft model. 

(3) Subroutine SUMMRY was dummied-out. This subroutine prints a very 
brief (one line per case) summary at the end of each computer run, 
which can be effectively deleted without degrading any program 
capabilities. This deletion eliminates some significant storage 
arrays. 

(4) To obtain the capability of receiving full program output in the 
event of execution time-out, the UNIVAC 1100 sponsored subroutine 

T I ML FT was added to the program. Also added was subroutine REMTIM, 
which acts as an interface between the HILTOP code and TIMLFT. 

The associated program input variable is ITF, which has a default 
value of 3 (seconds) pertaining to the IBM 360/91, and it may be 
necessary to increase this value (e.g., ITF = 10) on the UNIVAC 1100. 
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(5) Overlay segment LIMB0 was created for three subroutines (PARINC, 
PUNCH, and VPRINT) associated with program options which are 
rarely exercised. These options, which are then effectively dummied- 
out on the UNIVAC 1100, are described in the HILTOP program report 
by the Namelist inputs (respectively) KPART, MPUNCH and ALTITU. 
Suppression of these options does not degrade the program's basic 
electric propulsion mission analysis capability. Attempt to use 
any of these options on the UNIVAC could cause an abrupt execution 
halt due to an invalid overlay occurrence. 

The card images representing the HILTOP overlay tree on the UNIVAC 
1100 are displayed on the next page. 


HILTOP Overlay Tree for Uni vac 1100 


SEG R00T 

IN MFILE.A, .C, .DUMC0R, .L, .0,.Q,.R, .S s .U».AK,.AL,.AN 
IN MFILE.AQ..AR..AS, .AX,.REMTIM,.BB,.BO, .BE..DUMSUM, .BM 
IN MFILE .BP, .TIMLFT, .BQ, . VD0T, .VMAG, .BW 
SEG SEGQPR*, (R00T) 

IN MFILE.AV 

SEG INITAL*, (SEGQPR) 

IN MFILE.AH, .NEWINT. .AW, .BC, .BT 
SEG FINAL** (SEGQPR) 

IN MFILE.P, .DUM0RE, .AP 
SEG BALLIS** (SEGQPR) 

IN MFILE.AC, .AE,.AF 
SEG TRAJEC*, (R00T) 

IN MFILE.ALBED0, .E, .F,.G,.H,.J,.M, .N,.T,.V,.W,.X,.Y,.Z 
IN MFILE.AA,.AB,.AD,.AG,.AI,.AJ,.AT,.AY,.AZ,.BA,.BF,.BG,.BH 
IN MFILE.BK, .B0, .BR, .BS, .BX 
SEG SBEGIN*, (R00T) 

IN MFILE. D 

SEG LIMB0*, ( R00T ) 

IN MFILE. A0, . AU, ,BV 
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APPENDIX B 


h -h , 
o rof 


Outer Loop 


THRUST - false 
>f - large negative 


n - n + Aa 


Increment 

Number of thrusters 


Set optimal number 
of thrusts re and 
optimal voltage: 


V_ * Voltage {t ) 
I ^ ' save' 


Print message 
’INSUFFICIENT 
POWER TO THRUST” 


Inner Loop 


l- t + 1 


Increment 
Voltage Index 


Select Vj - Voltage (l). 
Compute beam current 1^ 

which satisfies the power 
constraint corresponding to 


n and V^ . 


I 

. max B 


Set! >0 

B max 


THRUST 


^orc&T\^ no 

s Tta*3tln^P ) 


h -h 
. o ref. 


Print message 
* 'REQUESTED FORCED 
THRUSTING NOT 
POSSIBLE rr 




THRUST = true 


THRUST 

\ / 


Compute Hamiltonian h 


Compute beam current I 

B 

which satisfies the power 
constraint corresponding to 


h -h 
ref 

n =• n 

save 

t = l 

save 


Print message 
,T ERROR ,r 




Set logical Indicators enabling exterlc*- logic to monitor 
the necessary Hamiltonian switch functions. 


Selection of Optimal Voltage and Number of Thrusters nt any Point In Time 
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